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Abstract 

The event-by-event fluctuations of suitably chosen observables in heavy ion colli- 
sions at SPS, RHIC and LHC can tell us about the thermodynamic properties of the 
hadronic system at freeze-out. By studying these fluctuations as a function of varying 
control parameters, it is possible to learn much about the phase diagram of QCD. As 
a timely example, we stress the methods by which present experiments at the CERN 
SPS can locate the second-order critical endpoint of the first-order transition between 
quark-gluon plasma and hadron matter. Those event-by-event signatures which are 
characteristic of freeze-out in the vicinity of the critical point will exhibit nonmono- 
tonic dependence on control parameters. We focus on observables constructed from 
the multiplicity and transverse momenta of charged pions. We first consider how the 
event-by-event fluctuations of such observables are affected by Bose-Einstein correla- 
tions, by resonances which decay after freeze-out and by fluctuations in the transverse 
flow velocity. We compare our thermodynamic predictions for such noncritical event- 
by-event fluctuations with NA49 data, finding broad agreement. We then focus on 
effects due to thermal contact between the observed pions and a heat bath with a 
given (possibly singular) specific heat, and due to the direct coupling between the crit- 
ical fluctuations of the sigma field and the observed pions. We also discuss the effect of 
the pions produced in the decay of sigma particles just above threshold after freeze-out 
on the inclusive pion spectrum and on multiplicity fluctuations. We estimate the size 
of these nonmonotonic effects which appear near the critical point, including restric- 
tions imposed by finite size and finite time, and conclude that they should be easily 
observable. 
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1 Introduction and Outline 



The goal of this paper is to motivate a program of heavy ion collision experiments aimed at 
discovering an important qualitative feature of the QCD phase diagram, namely the critical 
endpoint at which a line of first order phase transitions separating quark-gluon plasma from 
hadronic matter comes to an end HJ. The possible existence of such an end point, denoted E, 
in the temperature (T) vs. baryon chemical potential (//) plane has recently been emphasized 
and its universal critical properties have been described [El Ej|. The point E can be thought 
of as a descendant of a tricritical point in the phase diagram for 2-flavor QCD with massless 
quarks. As pointed out in [[[]], observation of the signatures of freezeout near E would 
confirm that heavy ion collisions are probing above the chiral transition region in the phase 
diagram. Furthermore, we would learn much about the qualitative landscape of the QCD 
phase diagram. 

In a previous letter [p]], we have laid out the basic ideas for observing the critical end- 
point. The signatures proposed in JTJ are based on the fact that such a point is a genuine 
thermodynamic singularity at which susceptibilities diverge and the order parameter fluc- 
tuates on long wavelengths. The resulting signatures all share one common property: they 
are nonmonotonic as a function of an experimentally varied parameter such as the collision 
energy, centrality, rapidity or ion size. The goal of the present paper is to develop a set 
of tools which will allow heavy ion collision experiments to discover the critical endpoint 
through the analysis of the variation of event-by-event fluctuations as control parameters 
are varied. 

Once experimentalists vary a control parameter which causes the freeze-out point in the 
(T, /i) plane to move toward, through, and then past the vicinity of the endpoint E, they 
should see all the signatures we describe first strengthen, reach a maximum, and then de- 
crease, as a nonmonotonic function of the control parameter. It is important to have a control 
parameter whose variation changes the fi at which the system crosses the transition region 
and freezes out. The collision energy is an obvious choice, since it is known experimentally 
that varying the collision energy has a large effect on \i at freeze-out. Other possibilities 
should also be explored.^ 

An example of nonmonotonic signatures in a different but analogous context is the rise 
and fall in the number of large fragments as a function of total observed multiplicity in 
multifragmentation experiments || in low energy nuclear collisions. These experiments allow 
us to confirm the existence and study the properties of another critical point — the end point 
of the first-order nuclear liquid-gas transition (boiling of the nuclear matter liquid to yield 
a gas of nucleons) || [5]. This point occurs at a temperature of order 10 MeV, much lower 
than the one we are studying |TJ. 

The analogy which is perhaps most familiar is with the phenomenon of critical opalescence 
observed in most liquids, including water. As the fluid cools down under conditions such that 
it passes near the end point of the boiling transition, it goes from transparent to opalescent to 
transparent as the end point is approached and then passed. This nonmonotonic phenomenon 

1 If the system crosses the transition region near E, but only freezes out at a much lower temperature, the 
event-by-event fluctuations will not reflect the thermodynamics near E. In this case, one can push freeze-out 
to earlier times and thus closer to E by using smaller ions.jl], [| 
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is due to the scattering of light on critical long wavelength density fluctuations, and thus 
signals the universal physics unique to the vicinity of the critical point. 

The universal property of systems in the vicinity of a second order critical point is the 
anomalous increase of thermodynamic fluctuations of the order parameter and related ob- 
servables. Here we consider a specific system, namely the hadronic matter created in a 
heavy ion collision at the time interactions freeze out. Our generic expectation is that the 
event-by- event fluctuations of suitable observables increase in the vicinity of a critical end- 
point. In this paper we calculate the magnitude of the resulting nontrivial effects, and make 
predictions which, we hope, will allow experiments to find the endpoint E. 

It is clear that before we can achieve this goal we must develop sufficient understanding 
of non- critical event-by-event fluctuations. Large acceptance detectors, such as NA49 and 
WA98 at CERN, have made it possible to measure important average quantities in single 
heavy ion collision events. For example, instead of analyzing the distribution of charged 
particle transverse momenta obtained by averaging over particles from many events, we can 
now study the event-by-event variation of the mean transverse momentum of the charged 
pions in a single event. The event-by-event variation of particle abundance ratios and even of 
the HBT radii are also becoming available. Although much of this data still has preliminary 
status, with more statistics and more detailed analysis yet to come, some general features 
have already been demonstrated. In particular, the event-by-event distributions of these 
observables are as perfect Gaussians as the data statistics allow, and the fluctuations — the 
width of the Gaussians — are small. [0 

This is very different from what one observes in pp collisions, in which fluctuations are 
large. These large non-Gaussian fluctuations clearly reflect non-trivial quantum fluctuations, 
all the way from the nucleon wave function to that of the secondary hadrons, and are not 
yet sufficiently well understood. As discussed in || ||, thermal equilibration in AA collisions 
drives the variance of the event-by-event fluctuations down, close to the value determined 
by the variance of the inclusive one-particle distribution divided by the square root of the 
multiplicity. In pp physics one can hope to extract quantum mechanical information about 
the initial state from event-by-event fluctuations of the final state; in heavy ion collisions 
equilibration renders this an impossible goal. In AA collisions, then, the new goal is to 
use the much smaller, Gaussian event-by-event fluctuations of the final state to learn about 
thermodynamic properties at freeze-out. 

What can we learn from the magnitude of these small fluctuations and their dependence 
on the parameters of the collision? Do they contain any more information than the corre- 
sponding moments of one-particle inclusive distributions? Some of these questions have been 
addressed in JTD, where it was pointed out that, for example, temperature fluctuations 



are related to heat capacity via 

((AT) 2 ) _ 1 m 
T 2 CV(T)' 1 1 

and so can tell us about thermodynamic properties of the matter at freeze-out. Similar ideas 



m 



11|] relate fluctuations in the occupation of certain momentum bins with d^i/dN and the 
average quantum density in phase space. Furthermore, Mrowczyhski has discussed the study 
of the compressibility of hadronic matter at freeze-out via the event-by-event fluctuations of 
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the particle number [12] and Gazdzicki[13] and Mr6wczyhski[14] have considered event-by- 
event fluctuations of the kaon to pion ratio as measured by NA49@. 

In this paper, we focus on observables constructed from the multiplicity and the momenta 
of the charged particles in the final state, as measured by NA49. It should be possible to use 
similar methods to analyze the event-by-event fluctuation of other classes of observables. For 
example, if it were possible to measure the baryon to pion ratio, analyses analogous to those 
we discuss would lead to the thermodynamic susceptibility d 2 fl/dfj 1 2 . As the neutrons are not 
observed, this analysis is not available. However, event-by-event fluctuations of the kaon to 
pion ratio may yield similar information. Another example is the data obtained by WA98 on 



the event-by-event fluctuation of the charged particle to photon ratio. JT^] They find a Gaus- 



sian distribution, and therefore constrain non-equilibrium processes in which long wavelength 
disorientations of the chiral condensate are excited, as these introduce non-Gaussianity. We 
leave the extension of the methods of this paper to the study of thermodynamic implications 
of the NA49 Gaussian distribution of event-by-event Kj n ratios and of the WA98 Gaussian 
distribution of event-by-event n /^ ratios for future work. 

Thermodynamic relations like ([I]) suggest the following strategy. Measure the mean 
transverse momentum of the charged pions in each event in an ensemble. Since the inclusive 
average of the transverse momentum of pions from an ensemble of events reflects (although 
does not equal) the temperature of the ensemble, perhaps one can use pt, the mean transverse 
momentum of the pions in a single event ,| as a proxy for the temperature of a single event, 
and so use (|1|) to obtain Cy. One of the lessons of the results we present below is that this 
strategy is too naive. To see a sign of this, consider another fundamental thermodynamic 
relation, namely that the event-by-event fluctuations of the energy E of a part of a finite 
system in thermal equilibrium are given by 

((AE) 2 )=T 2 C V (T) . (2) 

For a system in equilibrium, the mean values of T and E are directly related by an equation 
of state E(T); their fluctuations, however, have quite different behavior as a function of 
Cy, and therefore behave differently when Cy diverges at a critical point. So, is the Cy- 
dependence of the event-by-event fluctuations of pt like that of AT in ([[]) or like that of AE 
in (|2|)? We will show that pr fluctuations are not like either, although their Cy-dependence 
is more similar to that of AE, in the sense that the fluctuations of pr grow at the critical 
point. 

Most of our analysis is applied to the fluctuations of the observables characterizing the 
multiplicity and momenta of the charged pions in the final state of a heavy ion collision. 
There are several reasons why the pion observables are most sensitive to the critical fluctua- 
tions. First, the pions are the most numerous hadrons produced and observed in relativistic 
heavy ion collisions. A second, very important reason, is that pions couple strongly to the 
fluctuations of the sigma field (the magnitude of the chiral condensate) which is the order 
parameter of the phase transition. Indeed, the pions are the quantized oscillations of the 
phase of the chiral condensate and so it is not surprising that at the critical end point, where 



2 We denote the mean transverse momentum of all the pions in a single event by px rather than (pr) 
because we choose to reserve (. . .) for averaging over an ensemble of events. 
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the magnitude of the condensate is fluctuating wildly, signatures are imprinted on the pions. 
By Section 5, we will have built up the technology needed to analyze these signatures. 

Before we outline the structure of the paper, the following comment is in order. We 
assume throughout that freeze-out occurs from an equilibrated hadronic system. If freeze- 
out occurs "to the left" (lower /i; higher collision energy) of the critical end point E, it occurs 
after the matter has traversed the crossover region in the phase diagram. If it occurs "to the 
right" of E, it occurs after the matter has traversed the first order phase transition. This is 
the situation in which our assumption of freeze-out from an equilibrated system is most open 
to question. First, one may imagine hadronization directly from the mixed phase, without 
time for the hadrons to rescatter. Hadronic elastic scattering cross-sections are large enough 
that this is unlikely. Second, one may worry that the matter is inhomogeneous after the 
first order transition, and has not had time to re-equilibrate. Fortunately, our assumption 
is testable. If the matter were inhomogeneous at freeze-out, one can expect non-Gaussian 
fluctuations in various observables |16[] which would be seen in the same experiments that seek 
the signatures we describe. We focus on the Gaussian thermal fluctuations of an equilibrated 
system, and study the nonmonotonic changes in these fluctuations associated with moving 
the freeze-out point toward and then past the critical point, for example from left to right 
as the collision energy is reduced. 

Although our central point is the analysis of the critical fluctuations in the vicinity of the 
point E, we must first present an extensive analysis of the noncritical fluctuations, which are 
the background on top of which critical effects must be sought. Thus, this paper is organized 
as follows: Sections 2 and 3 analyze the background noncritical fluctuations; Section 4 
analyzes a particular (negative) contribution to the noncritical fluctuations which disappears 
near the critical point; Sections 5 and 6 analyze the critical fluctuations themselves. 

We begin in Section 2 by discussing the simplest case we can imagine, namely the fluc- 
tuations in an ideal Bose gas of pions. This allows us to establish some notation and to 
explain several conceptual issues. In particular, we explore the relation in this simplest 
case between the ensemble (i.e., event-by-event) variance and the variance of the inclusive 
one-particle distribution obtained by averaging over particles from many events. We also 
point out that the correlation between the multiplicity and an intensive observable, like the 
mean transverse momentum, only receives contributions from nontrivial effects such as Bose 
enhancement, energy conservation or interactions. This correlation is in general small, but 
we see in Sections 4 and 5 that it can increase by a large factor near the critical point. We 
derive results in Section 2 and throughout which are valid in the thermodynamic limit. In 
an Appendix, we explain the subtleties of constructing estimators for the relevant quantities 
using a finite sample of events each with a finite number of pions. 

Our goal in Section 3 is the inclusion of various effects neglected in Section 2, except that 
we continue to assume that freeze-out is not occurring in the vicinity of the critical point. 
We model the matter in a relativistic heavy ion collision at freeze-out as a resonance gas in 
thermal equilibrium, and begin by calculating the variance of the event-by-event fluctuations 
of total multiplicity. Our result suggests that about 75% of the fluctuations seen in the data 
are thermodynamic in origin. Our prediction is strongly dependent on the presence of the 
resonances; had we not included them, our prediction would have been significantly lower, 
farther below the data. 
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Fluctuations in extensive observables like the total multiplicity N are sensitive to non- 
thermodynamic variation in the initial size of the system which later thermalizes. Sources of 
such variation include: (i) the distribution of impact parameters; (ii) fluctuation in the initial 
positions of the nucleons; (iii) quantum fluctuations of the NN cross section jLjJ described by 
the wave function of the nucleon, which can be thought of as fluctuations in the effective size 
of the nucleons at the initial moment of the collision. All these effects lead to fluctuation in 
the number of spectator nucleons, and thus in the initial size of the interacting system which 
later thermalizes. We plan to evaluate the size of these contributions to fluctuations in N 
elsewhere. In this paper, we constrain the magnitude of these nonthermodynamic effects by 
comparing thermodynamic predictions for the fluctuations in N to the data. 

We then turn to a calculation of the variance of the event-by-event fluctuations of the 
mean transverse momentum, px- This is an intensive variable and should, therefore, be less 
sensitive to nonthermodynamic variations in the initial size of the system. We calculate 
numerically the thermodynamic contribution from "direct pions" , already present at freeze- 
out, and from the pions generated later by resonance decay. We include Bose effects and 
the effects of flow and find both to be small. We compare our results to those found by 
NA49 for central Pb-Pb collisions at 160 AGeV, and find broad agreement. We do not 
attempt to include purely experimental effects, such as those due to two-track resolution, 
and so do not expect precise agreement. Our goal is to compare observed variance with 
thermodynamic expectations and to see whether they are consistent. Our results support 
the general idea that the small fluctuations observed in AA collisions, relative to those in pp, 
are consistent with the hypothesis that the matter in the AA collisions achieves approximate 
local thermal equilibrium in the form of a resonance gas. Once data is available for other 
collision energies, centralities or ion sizes, the present NA49 data and the calculations of this 
section will provide an experimental and a theoretical baseline for the study of variation as 
a function of control parameters. 

In Sections 4 and 5, we analyze how the proximity of the critical endpoint to the freeze- 
out point is reflected in the fluctuations. We begin in Section 4 by making the idealization 
that the pions which one observes are an ideal Bose gas in thermal contact with a heat 
bath which includes the sigma field. The heat capacity of this heat bath is therefore infinite 
at the critical point. This treatment neglects the aim coupling, which allows the critical 
fluctuations of the sigma field to influence the pion fluctuations directly, rather than just by 
thermal contact. 

The dominant effects of the critical fluctuations on the pions are the direct effects oc- 
curring via the ami coupling. The idealization of Section 4 is nevertheless useful, because 
it allows us to explain and illustrate an important point not made clear in M related to 
the practical application of ([!]) • The fluctuations of the temperature depend on what "me- 
chanical" observable (such as the energy, for example) is measured, and how the measured 
observable is converted into a temperature. In particular, these fluctuations depend on what 
part of a system is used as a thermometer. Eq. (H) describes a particular case when the whole 
system of interest is used as a thermometer. It requires us to use the equation of state, T(E), 
of the whole system of interest to translate the energy, which is measured in this case, into 
the temperature [ IS| . The fluctuations of "mechanical" variables, such as energy, increase at 
the critical point, as in (0). Because T(E) is singular at the critical point, the fluctuations 
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of T decrease, and vanish at the critical point where Cy — > oo. It is a fact that what we 
measure are the mechanical observables, and since we in general only know T(E) for simple 
systems we call thermometers, we cannot apply (|1]) to the complicated system of interest. 
We illustrate these points by evaluating the fluctuations of several observables in an ideal gas 
of detected pions (the thermometer) which is in thermal contact with an undetected non- 
ideal, possibly singular, heat bath. The effect we find vanishes at the critical point, where 
the specific heat of the heat bath diverges due to the fluctuations of the sigmas therein, and 
so provides a nonmonotonic signature. The effect involves a reduction in the fluctuations 
of the mean transverse momentum of the pions. What makes it distinctive is that it also 
involves an anti-correlation between fluctuations of pion occupation numbers with different 
momenta. We show that this phenomenon follows directly from energy conservation, and 
conclude that it is much more robust than the idealizations we use to describe it. This 
signature is present when the system freezes out far from the critical point, and is reduced 
near the critical point. It should be observable in present data on central PbPb collisions at 
160 AGeV, even if freeze-out is not occurring near the critical point in these collisions. 

Section 5 describes what we believe to be the dominant event-by-event signatures directly 
related to the divergent correlation length which characterizes the critical point. We apply 
much of the technology built up over the preceding sections in Section 5 to study the effect 
of the interaction of the pions with the almost classical sigma field. We find a large increase 
in the fluctuations of both the multiplicity and the mean transverse momentum of the pions. 
This increase would be divergent in the infinite volume limit precisely at the critical point. We 
apply finite size and finite time scaling to estimate how close the system created in a heavy ion 
collision can come to the critical singularity, and consequently how large an effect can be seen 
in the event-by-event fluctuations of the pions. We conclude that the nonmonotonic changes 
in the variance of the event-by-event fluctuation of the pion multiplicity and momenta which 
are induced by the universal physics characterizing the critical point can easily be between 
one and two orders of magnitude greater than the statistical errors in the present data. 

Once we have analyzed the effects of the sigma field on the fluctuations of the pions, in 
Section 6 we ask what becomes of the sigmas themselves. Assuming that freeze-out occurs 
near the critical point, they are numerous at freeze-out and they can only decay later, once the 
sigma mass has risen above twice the pion mass. This results in a nonmonotonic signature 
of the critical point which can be observed even without an event-by-event analysis. We 
calculate the momentum distribution of these low momentum pions produced in the delayed 
decays of the sigmas. We close by analyzing the enhanced event-by-event fluctuations of the 
multiplicity of these low momentum pions. 

We end the paper with a summary of the different contributions to the event-by-event 
fluctuations which we have analyzed, and a more general look to the future. In striving 
to provide analyses which will assist experimentalists to use the universal properties of the 
critical point to learn its location, we hope that we have in addition provided a set of tools 
for event-by-event analyses of heavy ion collisions which will prove useful in the study of the 
thermodynamics of QCD in a variety of contexts in the future. 
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2 Thermodynamic Fluctuations in an Ideal Bose Gas 



2.1 The Basics 

We begin by recalling text-book facts about the thermodynamics of an ideal Bose gas which 
are relevant to our event-by-event analysis. Little in this section is new, but it is nevertheless 
a very helpful exercise and will allow us to establish some notation. The basic fact is that 
every quantum state of a system of identical spinless Bose particles is completely character- 
ized by a set of numbers, n p — the occupation numbers for the one-particle states labeled 
by momenta p. All thermodynamic quantities are functions of these numbers and thus all 
we need to know is the fluctuations of n p from one member of the ensemble (one event) 
to another. The first step toward a characterization of these fluctuations is the ensemble 
average of the occupation number for the mode with momentum p, namely 

^ = e e P /T _ V ( 3 ) 



where e p = u p — \i and, as usual, u p = \Jp 2 + m 2 . Next, we need the deviation, An p = 
n p — (n p ) , whose mean square average in the ensemble is given by 

{{Anp)2) = T ^= k - 1) 2 = + K)) " * (4) 

We have introduced notation v p for this quantity which will be used frequently below. This 
expression is "microscopic" , in the sense that it is written for a single mode in momentum 
space. However, it can be derived "macroscopically" as follows. The fluctuations in the total 
particle number 

N = Y,n p (5) 

p 

are given by[[H| 

«A*)') = r(f) T . (6) 

Because the fluctuations of different modes are statistically independent, we can elevate this 
relation to the microscopic form (f|), and indeed to 

(An p An k ) = ((An p ) 2 )S pk = v 2 p 5 pk . (7) 

The correlator (An p Arik) is the central quantity which we will calculate repeatedly through- 
out this paper, as we proceed beyond the ideal Bose gas. 

The correlator in (0) enters in the calculation of the event-by-event mean square deviation 
of any generic thermodynamic variable of the form 

Q = S % u p ■ ( 8 ) 
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Indeed, since AQ = Q — (Q) = J2 P QpAn p , we find that 

((AQ) 2 ) = £ g p g fc (An p An fe > = £ g> 2 . (9) 

pk P 

The quantity Q could be the total energy 

£ = £n p e p , ((AE) 2 ) = J244-i ( 10 ) 

p p 

or it could be the total transverse momentum, J2 P {PT)p n p] or simply the total particle number 

iV = 5> p , <(AiV) 2 >=5> 2 . (11) 

p p 

For future reference we also give here an expression for the heat capacity Cy of the Bose 
gas at constant V and \i: 

*- T (£).~ T (£l 

Using the expression for the thermodynamic potential: 



tt = Tj2H l ~ e~ ep/T ). (13) 
p 



^ -«»/T 

P 

one finds 

J^2 



Cv = ^E44- ( 14 ) 

p 



Comparing to (|T0| ) we find the well-known relation 

((AE) 2 ) = T 2 C V , (15) 
which is valid for any system in equilibrium. 

2.2 Energy per Particle: Event-by-Event Average vs. Single- 
Particle Inclusive Average 

Let us now compute the fluctuation of an intensive observable, such as the mean energy per 
particle e = E/N, where E and iV are extensive, or px, the mean transverse momentum 
per particle in a single event. Analyzing the member-of-the-ensemble-by-member-of-the- 
ensemble fluctuations of the mean energy per particle in a single member of the ensemble is 
a good warmup. We henceforth begin to refer to members of the ensemble as events. For 
small fluctuations (and AE/E ~ N~ l l 2 <C 1 is small) we can write 





E 


(AE 


AN\ 


(#)■ 




K E 


N ) 
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Ai'-H =^((AE) 2 +(^)\aN) 2 -2(^)aEAn). (17) 



Now, we square: 

Then we average. We already know ((AE) 2 ) and ((AN) 2 ), but we also need 

(AEAN)=Y,t P v 2 p , (18) 

p 

which is obtained in the same way as before. Putting this all together, we find 

((A e ) 2 ) = ^E(^-(^))% 2 - (19) 

Let us now compare (|T9|) to the variance of the inclusive single particle average energy 
per particle. To do this, we introduce the notation 

r nc = E / Efo> = EfcW/W • ( 2 °) 
p p p 

Whereas (...) denotes an average over events of some property of a single event, "ttt^ 1110 denotes 
an inclusive average of a property of a single pion over all pions in the ensemble of events, 
without reference to in which event each pion occurs. It is more convenient for theoretical 
purposes to work with occupation numbers n p , and the inclusive average is then done n p 
by n p as defined in (^). The subscript p on the left hand side of ( j2"0| ) reminds us that 
the average is done momentum bin by momentum bin: it is q p which is being averaged, 
not q = Q/N. However, the quantity % inc is p- independent. Were we only interested in a 
quantity like (e), there would be no need to take care with definitions because averaging a 
single particle quantity pion by pion is the same as first averaging it over an event, and then 
averaging event-by-event: 

(e) = (E/N) = ^ c . (21) 

This is not true for fluctuations about the mean, as we see by using our definitions to rewrite 
(0) as 

<(A e ) 2 > = ^y( ep -(6)) 2 (l + K)) mC . (22) 

The same formula holds if e = E/N is replaced by any quantity of the form q = Q/N, for 
example by the mean transverse momentum per event. 

The lesson we learn from fl22"|) is that up to the Bose enhancement factor (1 + (n p )), the 
ensemble (alias event-by-event) fluctuations of intensive quantities, such as the energy per 



— inc 

particle, are indeed given by the variance of the single particle distribution (e p — (e)) 2 and 
the central limit theorem which dictates the factor 1/(N). We see that the effect of the Bose 
factor is to increase the variance of the event-by-event distribution relative to that of the 
inclusive distribution. 

When we apply formulae like those we have just derived which are valid in the thermody- 
namic limit to heavy ion collision data, we will need to construct estimators for the relevant 
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quantities using a finite ensemble of events, in which the number of particles in each event 
is also finite. We describe how this should be done in an Appendix. 

Having discussed the fluctuations of extensive and intensive quantities, we end this Sec- 
tion by considering the cross correlation between an intensive observable and the extensive 
observable N. For example, let us calculate (AeAN). Using ingredients we have spelled out 
above, we find 

(A^AiV) = ^yX>J 0* - (c» = EW 2 (h - <e» . (23) 

Note that the terms proportional to (n p ) have cancelled, and the remaining term, propor- 
tional to (n p ) 2 , is obviously due to the Bose effect. This result applies to any such correlation; 
for example we could have used px instead of e. The lesson we learn is this: cross corre- 
lations between iV and intensive observables are generally small, because they receive no 



contribution if one takes the classical ideal gas limit. Recall that in (|21f) we find a dominant 
contribution to the event-by-event variation coming from the variation of the inclusive sin- 
gle particle distribution. In (Iffi), this effect cancels and the remaining effects due to Bose 
enhancement dominate. This means that the nontrivial effects on the pions due to their in- 
teractions and due to energy conservation and thermal contact with other degrees of freedom 
only need to be larger than the effects of Bose enhancement in order to dominate this cross 
correlation. 



3 Noncritical Thermodynamic Fluctuations 
in Heavy Ion Collisions 

In this section we proceed to quantitative estimates of the magnitude of noncritical event-by- 
event fluctuations in heavy ion collisions. As an example of an extensive quantity we use the 
total charged pion multiplicity of an event; as an example of an intensive quantity we use the 
mean transverse momentum pt of the charged pions in an event. We compare some of our 
estimates to preliminary data from the NA49 experiment at CERN SPS on PbPb collisions 
at 160 AGeV, and find broad agreement. In this section, and throughout this paper, we 
assume thermal equilibrium at freeze-out. In this section, but not throughout this paper, 
we assume that the system freezes out far from the critical point in the phase diagram, and 
can be approximated as an ideal resonance gas when it freezes out. The results obtained 
seem to support the hypotheses that most of the fluctuation observed in the data is indeed 
thermodynamic in origin and that PbPb collisions at 160 AGeV do not freeze out near the 
critical point. 

3.1 Pion Gas at Thermal Freeze-out and Bose Enhancement 

The observed spectrum of pions reflects the distribution of pion momenta at the time of ther- 
mal freeze-out, namely the time at which the interaction rates fall behind the expansion rate. 
After this time, one can approximately neglect energy /momentum exchange interactions and 
consider the momenta of particles as frozen. Freeze-out is by definition the time at which the 
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system ceases to be in thermal equilibrium. However, if the system has thermalized before 
it freezes out, then even after freeze-out one has a thermal distribution of pion momenta,^ 
approximately with a single temperature over the whole system. This standard idealization 
at this point seems sufficient to describe the data. (Particles which interact more weakly 
than pions freeze out earlier, at a higher temperature. We leave such particles, together with 
details of the dynamics of the freeze-out of the pions, to future work.) 

We start with the simplest model for the pions at freeze-out — the ideal Bose gas. This 
allows us to use the results of the previous section. Later in this section, we add pions 
produced by the decay of resonances as well. The isospin degeneracy of the pions requires 
a small modification to the formulae of the previous section. Since only the momenta of 
charged pions are observed, we must only count 7r + and 7r~. Because 7r + and 7r~ are distinct, 
the Bose enhancement factor is reduced from 1 + n p to 1 + (n p /2), where n p counts the total 
number of charged pions. This is the consequence of the fact that only identical pions can 
interfere|] 

We begin by showing that the effect of the Bose enhancement is very sensitive to a nonzero 
pion chemical potential fi n (not to be confused with the baryon number chemical potential 
fi). Let us first remind the reader why a nonzero fi n may be needed. The pion chemical 
potential is not a thermodynamic conjugate to any fundamentally conserved quantity, and 
is the same for pions of all charges. It is supposed to represent the over-population of pion 
phase space. It allows for the possibility that even though the momenta of the pions are 
in equilibrium at freeze-out, their number density is not. This arises because all reactions 
which can change the number of particles, and thus keep this quantity in equilibrium, have 
small cross sections at the relevant low energies. In contrast, elastic re-scattering is strongly 
enhanced by resonances (such as A,N* for ttN, a, p for tttt etc). As a result, thermal 
equilibrium of momenta is maintained to a lower freeze-out temperature, whereas chemical 
freeze-out (below which particle numbers do not change) occurs somewhat earlier. There 
is therefore a window of time between chemical and thermal freeze-out during which the 
system evolves with fixed pion number; during this time a pion chemical potential naturally 
develops. At chemical freeze-out, fi n = 0. As the temperature then continues to drop while 
the pion number remains fixed, p n increases. For an overview of pion kinetics and references 
see Practical calculations of the magnitude of the effect for heavy ion collisions at 

CERN SPS can be found in pD| . The conclusion inferred from this analysis is that the pions 
in central PbPb collisions at SPS energies freeze out at a temperature Tf « 120 MeV with 
H„ « 60 MeV. 

Now we return to the calculation of the Bose enhancement of fluctuations of some generic 
single-particle intensive observable q = Q/N. If we use the notation v e b e for the event-by- 
event variance and v- mc for the variance of the inclusive distribution: 



<be(?) = <(Ag) 2 > and < c (g) = (q p - , (24) 



3 As is very accurately the case for the cosmic microwave background radiation, ten billion years after its 
freeze-out. 

4 It is easy to see that (An^An^) = (n p )(l + (n l p ))8 l: > 5 p k, where i,j = +, — . On the other hand, from 
n p = J2i n P ^ follows that: (n p ) = (n p )/2, and (An p An k ) = (n p )(l + (n p )/2)S p k. 
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and define the ratio 



F _ (N)vlM 
then we can write the result (^) for the ideal Bose gas as 



B 



1 + 



l Jd 3 p(q P -(Q)) 2 {n P ) 2 
2 Jd3p(q p -(q)y(n p ) 



(25) 



(26) 



The factor of 1/2 appears because, as discussed above, there are two charged pions. As 
we consider effects not present in an ideal gas, we will find that the ratio F is not given 
simply by the Bose enhancement factor Fb- It is a product of Fb and other factors which 

we estimate later in this section and in subsequent sections. 

The dependence of Fb on is shown in Fig. 1 for q = px = ^Jpl + Py. (Note that vf nc 
does depend on the pion chemical potential as well.) We have also shown Fb for pions with 
Pt < 300 MeV to demonstrate that restricting the acceptance to low-energy pions results 
in a larger Bose enhancement effect. It is worth noting that for more central collisions the 
thermal freeze-out temperature, Tf, is lower because the system is larger and freezes out 
later [|]]. Therefore, fi n should be somewhat larger and the Bose enhancement effect should 
somewhat increase event-by-event fluctuations for more central collisions. To conclude, the 
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Figure 1: Bose enhancement factor Fb = (N)v^ he (pT) / vf nc (px) describing the contribution of 
the Bose effects to the fluctuations of the mean transverse momentum in an ideal Bose gas 
of pions. Fb is plotted as a function of the pion chemical potential in MeV. The dashed 
line shows the Bose enhancement factor if only pions with low momentum px < 300MeV are 
included. 



Bose enhancement effect is sensitive to fi^, and leads to an increase in t> e b e by a factor of 
y/Fs, which typically results in an increase of the order of a few percent. 
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The effects of Bose enhancement on the variance of the fluctuations of pt in an ideal 
Bose gas have been considered previously^!], [7]]: our results are in quantitative agreement 



with theirs. These authors use the quantity 

% T = {N) l ' 2 v ehe (p T ) - v inc (p T ) (27) 

introduced in Ref. § as a measure of the Bose enhancement effect. As we discuss in 
the Appendix and below, when (N) is finite one must take care in defining $ PT to use an 
appropriate estimator for t> e be(Pr)- To compare our results with theirs, note that what we 
describe as y/Fs — 1-01 corresponds to $ PT = (0.01)t> inc (pr)- This already hints at what 
we will see below, namely that whereas $ PT depends on the flow velocity through fi nc (pr), 
ratios like Fb are much less sensitive to the effects of flow and are therefore more easy to 
calculate. 



3.2 Contribution of Resonances 

The hadronic matter produced in a heavy ion collision is not simply an ideal gas of pions. A 
number of approaches to heavy ion collisions have successfully treated the matter at freezeout 
as a resonance gas in thermal equilibrium. Our analysis of the fluctuations observed in 
present data lends support to this idea. 

Even the global properties of hadronic matter are strongly affected by resonances. Al- 
though at temperatures of interest (T < T c ) the Boltzmann factor exp(— M/T) for each 
resonance is small, it is partially compensated by the pre-exponential factor due to the large 
number of states involved. One may recall here the Hagedorn conjecture, that at a certain 
temperature the contribution to the energy density due to the resonances would diverge be- 
cause of the exponential growth of the density of states. Although this does not happen in 
practice, because the chiral phase transition occurs at a much lower temperature than any 
putative Hagedorn transition, one nevertheless finds that when relevant resonances are 
included, the energy density and pressure increase rapidly with T, and can be fitted by a 
power law 

e(T) ~ P(T) ~ T K (28) 

with the power k ~ 6 at zero baryon density^ One also finds that the heat capacity, 
normalized to the number of pions (which means that we have "decayed" all resonances, 
counting each p meson as 2 pions, each a; as 3, etc.) is 

{Cy / ^-k) resonance gas ~ 23, (^9) 

at Tf = 120 MeV, while for the ideal Bose gas of pions one has only: 

(CV /Anions « 14 (30) 

at the same temperature. This general observation [|ll|] already suggests that the resonances 
may affect the fluctuations considerably. 



5 For nonzero baryon density this effective power is even larger, but we will ignore this since it is only 
important at much higher baryon densities (and lower collision energies) than achieved at the CERN SPS. 
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The resonances play another role in the problem. Those which are present at freeze-out 
decay after freeze-out, and by definition this means that the pions they produce cannot 
rescatter. The pions observed in the data are therefore a sum of (i) "direct pions" which 
were pions at freeze-out and (ii) pions produced from the decay of resonances after freeze- 
out. Note that the direct pions also originate from resonances, in the sense that most of 
the low energy rescattering which occurs before freeze-out occur via resonances. As we 
have discussed above, inelastic reactions which change the number of participants freeze-out 
earlier than elastic scattering. What this means is that the multiplicities of the pions and 
resonances, although thermal, should be fixed not at the thermal freeze-out (Tf ~ 120 MeV 
in PbPb) but earlier, at chemical freeze-out. Fits of SPS data on ratios of particle yields to 
thermal models yield T ch ~ 160 - 170 MeV fl23jf| 

In the remainder of this subsection, we investigate three effects of the resonances on 
the event-by-event fluctuations of the extensive observable N, the number of charged pions 
per event, and the intensive observable pt, the mean transverse momentum per event. We 
first describe all three effects briefly, and then describe the simulation which we have used 
in order to investigate them. The first (and largest) effect is a direct contribution to the 
fluctuations of N, and indeed to any extensive observable. The multiplicity fluctuations in 
a classical ideal gas are characterized by ((AN) 2 )/(N) = 1 and for an ideal pion gas this 

ratio is 1 + (l/2)(n p ) mc due to Bose effects and is therefore a few percent larger than 1. 
This ratio is significantly larger for a resonance gas. Each resonance decays into several 
pions (for example, p —>■ 2ir,uj,r] — > 3n, etc), and this means that when the resonances 
decay after freeze-out, they significantly modify the statistics of pion number fluctuations. If 
the resonances themselves are produced randomly]], with a Poisson multiplicity distribution, 
their decay products are not. For example, if there were no direct pions and only one species 
of resonance which always decayed into d charged pions, the pions produced in this ensemble 
would have ((AN) 2 )/(N) = d. 

Resonances also affect the fluctuations of intensive observables, like the energy per pion or 
mean transverse momentum p?- The second effect we analyze arises because pions produced 
in resonance decays have a single-particle momentum spectrum which is similar but not 
identical to the thermal spectrum for the direct pions. The products of resonance decay 
populate the lowpr region of the spectrum somewhat more.0 In order to estimate u e be(pr), 
we must therefore include the change in v inc (pT) introduced by the pions produced by the 
decay of resonances after freeze-out. 

The third effect of resonance decays is that they contribute additional kinematic correla- 
tions between their decay products, which then have no chance to rethermalize. New terms 
arise in the correlator (|7|) which describes the fluctuations microscopically. For example, 
two-body decay (such as p — > tt + ti^ at rest) generates a term in the correlation function: 

(An+An fc ) = 5 p> - k C(p), (31) 

6 Note that this number is close to the critical temperature obtained in lattice simulations with zero baryon 
density Q. 

7 All resonances are heavy enough that Bose enhancement for them can be neglected. 
8 This effect is qualitatively similar to the effect of a nonzero pion chemical potential. A clear distinction 
between these effects in data analysis is still lacking. 
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where C(p) is proportional to the square of the pion fraction originating from p decays. 
The result of all such terms will be a change in twO^t) which can be parameterized as a 
new contribution F res to the ratio F of (p5|). (That is, we now have F = FBF res .) Instead of 
attempting to study all contributions like (|3l|) one by one, we address this effect and the first 
two by doing a simulation. We will see that the second and third effects we have described 
are both small. 

We have simulated a gas of pions, nucleons and resonances in thermal equilibrium at 
freeze-out, including the tt, K, rj, p, u, rf, N, A, A, £ and H, and then simulated the 
subsequent decay of the resonances. That is, we have generated an ensemble of pions in 
three steps: (i) Thermal ratios of hadron multiplicities were calculated assuming equilibrium 
ratios at chemical freeze-out. Following p3| , the values T ch = 170 MeV and /ibaryon = 200 



MeV were used, (ii) Then, a program generates hadrons with multiplicities determined 
at chemical freezeout, but with thermal momenta as appropriate at the thermal freeze-out 
temperature, which we take to be Tf = 120 MeV, with p n = 60 MeV. The last step (iii) 
is to decay all the resonances, using the appropriate subroutine from RQMD.0 Under these 
conditions, more than half of the observed pions come from resonance decays. 

We evaluate the variance of the fluctuations of the multiplicity of the pions obtained from 
the resonance gas as follows. For each species in the resonance gas, we label the different 
decay modes by an index i, and refer to the branching ratios for the species r as b\. For each 
decay we define d l r , the number of charged pions produced. From the simulation, we obtain 
the multiplicity of each resonance, N r . The total multiplicity of pions is N n = J2r,idlb l r N r 
and the multiplicity fluctuations are described by 

(N w ) E r diKN r ' ' 1 ' 

Bose enhancement increases this to ((AN^) 2 ) / (N^) pa 1.5.^] We see that the resonances 
increase the multiplicity fluctuations by a large factor, relative to the fluctuations of the 
direct pions alone. We compare this result to what is seen in NA49 data below. 

We now turn to the resonance-induced contribution to the fluctuations of the intensive 
observable px, which is much smaller. We begin with the effect on v- mc . Table 1 describes 
the single-particle inclusive distribution obtained from the simulation, assuming uncorre- 
lated particles in an equilibrium resonance gas at freeze-out. We see that the resonances 
change v- inc (pT) / (pr) only by a few percent. The contributions of correlations induced by 
resonance decay and of Bose enhancement to F are not included. The effects of flow are 
not included. We now discuss each in turn, and find that all yield small contributions to 
(iV) 1 / 2 u e be(pT) / '{Pt) relative to v inc (pT) / (pt) which we have evaluated in the table. 

We have estimated F Tes by slicing up the pions from Table 1 into varying numbers (up to 
2500) of events, and evaluating F. Since Bose enhancement is not included in the simulation, 

9 We treat particles which decay by weak interactions as stable, which raises an additional issue. Exper- 
imentally, some weak decays happen so quickly that they feed up into the observed pion spectra. We treat 
these particles as stable here; we hope that experimentalists make the appropriate corrections to the data. 

10 Event-by-event fluctuation in the resonance multiplicities N r , as may be computed, for example, in 
dynamical models in which the resonances themselves are produced by decays of "clusters" , may result in a 
small further increase in this ratio. 
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no. of pions, 
10 3 


(PT>, 

MeV 


MeV 


Vmc(pT)/{pT) 


"direct pions" only 


541 


283 


189 


0.67 


pions from resonances only 


713 


271 


177 


0.65 


all pions 


1254 


276 


183 


0.66 



Table V. Results of a numerical simulation of a resonance gas. The results include the 
effects of the the correlations induced by resonance decays on the inclusive pt spectrum. 
The simulation itself does not include Bose enhancement effects, and so can be thought of 
as the simulation of a single event with 1.254 • 10 6 pions with a mean p? of 276 MeV, or can 
be sliced up into smaller events. 



the F so obtained is just F ves . We find no statistically significant contribution to F, and 
conclude that |F res — 1| < 0.01. 

We now use the results of Section 3.1 to incorporate Bose enhancement effects, after 
noting the connection between Bose enhancement and resonance decay pions. There can be 
quantum interference between direct pions and resonance decay pions, or among resonance 
decay pions. It is well known that all resonances can be approximately separated into two 
groups: those which are short-lived and those which are long-lived. The former (e.g., p and 
A) have lifetimes much shorter than the duration of pion radiation from the fireball (i.e., 
the time over which freeze-out occurs). Therefore, their decay products interfere with other 
pions. The decay products of long-lived particles (e.g., u and rj) can only interfere with other 
pions if one selects pions with a very small energy difference \uj\ — u)2\ « T, where V is the 
width of the resonance. This is essentially impossible, and pions produced in the decay of 
long-lived particles therefore do not contribute to the Bose enhancement factor. This means 
that Fb — 1 should be multiplied by a factor (/direct + f P + /a + ...) 2 where the /'s are the 
fractions of all 7r mesons coming from short-lived sources. This same fraction enters the HBT 
correlation function, and is about 0.5. ||25|| So, we take Fb — 1 = 0.073 for p n = 60 MeV from 



Fig.(^) and reduce it by a factor of 0.5 yielding 

F B = 1.037 , (33) 

and therefore conclude that the effect of Bose enhancement is a small increase in v e h e {pr) by 
a factor of ^JFb~ = 1-018. 



3.3 The Effects of Radial Flow 

To this point, we have calculated the fluctuations in px as if the matter in a heavy ion 
collision were at rest at freeze-out. This is not the case: by that stage the hadronic matter is 
undergoing a collective hydrodynamic expansion in the transverse direction, and this must 
be taken into account in order to compare our results with the data. A very important point 
here is that the fluctuations in pion multiplicity are not affected by flow, and our prediction 
for them is therefore unmodified. Fluctuations in multiplicity ratios (e.g. K/tt) would also 
be unaffected. However the event-by-event fluctuations of mean pt are certainly affected 
by flow. The fluctuations we have calculated pertain to the rest frame of the matter at 
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freeze-out, and we must now boost them. A detailed account of the resulting effects would 
require a complicated analysis. Here we shall use the simple approximation p6| that the 
effects of flow on the pion momenta can be treated as a Doppler blue shift of the spectrum: 
n ipr) — > n(pr\/l — + P)- This blue shift increases (pr), and increases v^ipr), but 

leaves the ratio i>i nc (px)/(pr) (and therefore the ratio f e be(pr)/(pr)) unaffected. This ratio 
(the fourth column in Table 1) is therefore a good quantity to compare to experimental data, 
since our goal here is to extract information about thermodynamics and not about flow. 

However, event-by-event fluctuations in the flow velocity (3 must still be taken into ac- 
count. This issue was discussed qualitatively already in [[□]], where it was argued that this 
effect must be relatively weak. Here we provide the first rough estimate of its magnitude. 
The magnitude of the flow velocity is proportional to the integral of the pressure gradi- 
ent over the expansion time. Thus the fluctuations of the flow velocity are determined by 
the pressure fluctuations. The size of the pressure fluctuations is related to the adiabatic 



compressibility by the standard thermodynamic relation [18 



« AP ) a > = - T (^) ■ ( 34 ) 

For the resonance gas equation of state this gives 
((AP) 2 ) k T k 2 1 



P 2 K-1PV K-1S' 



(35) 



The entropy per pion in the ideal gas is around 2.4, and is larger for the resonance gas. We 
shall take S ~ 3N n for our estimate. 

Using the "blue shift" approximation we can write 



PT^PTf^ (36) 

where p r ^ st is the corresponding momentum in the rest frame of the matter. The fluctuations 
of the observed momentum are then related to the fluctuations in the rest frame (calculated 
above) and the flow velocity fluctuations through 

«WUM + ((Am , (37) 

where we have neglected corrections which are suppressed by 0(/3) relative to ((A/3) 2 ). The 
fluctuations in the flow velocity are given by 

32 ((AP) 2 ) fr„ 



PF) = ^ — ■ (38) 

The last factor on the right-hand side appears because the final velocity is proportional 
to the time integral of the pressure gradient over the entire evolution prior to freeze-out, and 
this integral is a sum over uncorrelated fluctuations in time. In a resonance gas one can 
discuss the typical duration of a collision (the lifetime of a typical resonance), and the time 
between collisions (the inverse of the scattering rate). Both are close to the "microscopic" 
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time scale r m i cro ~ 1 fm/c. The expansion duration relevant for radial flow is actually 
much longer, r now »s 10 — 20 fm/c for central Pb-Pb collisions. This means that for each 
microscopic volume element one first does the time integral and obtains a "random walk" 
factor (rmicro/Vfiow) 1 / 2 ~ 1/4 in A/3. Then, the sum over uncorrected volume elements leads 
to a 1/ y/V or 1/ yN, which we have already seen in the expression ( pop for (AP) . 

The flow velocity can be estimated for our purposes from the ratio of (p^ st ) ~ 276 MeV 
given in Table 1 and the experimental (px) ~ 376 MeV observed by NA49. Thus, (3 ~ 0.3. 
Finally, putting all the estimates into eq. ( j38|) we find 

N((Af3f) « (0.1) 2 . (39) 

Note that although our estimate is uncertain at various points, the result is very small. 
Even if we have underestimated the size of ((A/3) 2 ) by a factor of 4, the contribution to 
VmciPr) I '{Pt) would only be 0.02. It is quite clear that the great bulk of v- mc J (p?) is thermo- 
dynamic, with the contributions of the fluctuations in the flow velocity being negligible in 
comparison. 

The largest uncertainty in our estimate for fine (pt)/(pt) is not due to the fluctuations in 
the flow velocity, which can clearly be neglected, but is due to the velocity itself. The blue 
shift approximation which we have used applies quantitatively only to pions with momenta 
greater than their mass[P6|. Because of the nonzero pion mass, boosting the pions does not 



actually scale the momentum spectrum by a momentum independent factor. Furthermore, 
in a real heavy ion collision there will be a position dependent profile of velocities, rather 
than a single velocity (3. A more complete calculation of v- mc (jpr)/ (pr) would require a better 
treatment of these effects in a hydrodynamic model; we leave this for the future. 

We obtain our final estimate of the magnitude of the event-by-event fluctuations of 
the intensive quantity px far from the critical point as follows. Using the estimate of 
v\nc(j>T) rest I '(pr) rest from Table 1 and Eqs. (^,^), we estimate that fluctuations in the flow 
velocity increase v- mc {pT)/ (pt) from 0.66 to 0.67. Multiplication by \/F~b then yields 

V 7 . ^ J « 0.68 , 40 
(Pt) 

subject to the uncertainties introduced by the blue shift approximation. 



3.4 Comparison with NA49 Data and Outlook 

In this section we compare our results with the NA49 data from central Pb-Pb collisions J7| 
summarized in Table §. As a first qualitative check of the predictions of our resonance 
gas model, we can look at the multiplicity fluctuations. It is clear that with no cut on 
centrality, one would see a very wide non-Gaussian distribution of multiplicity determined 
by the geometric probability of different impact parameters b. Gaussian thermodynamic 
fluctuations can only be seen if a tight enough cut in centrality is applied. The event-by- 
event iV-distribution found by NA49 when they use only the 5% most central of all events, 
with centrality measured using a zero degree calorimeter, is Gaussian to within about 5%. 
This cut corresponds to keeping collisions with impact parameters b < 3.5 fm.0] The non- 
Gaussianity could be further reduced by tightening the centrality cut further. We now ask 
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number of events 


98426 


total number of charged particles 


26587685 


(N) 


270.13±0.07 


Vebe(N) 


23.29 ± 0.05 


(Pt) 


376.75 ± 0.06 MeV 




282.16 ± 0.04 MeV 


Vebe(PT) 


17.27 ± 0.03 MeV 



Table 2: Preliminary NA49 dataj7|. The charged particles are taken from the kinematic 
region 0.005 < px < 2 GeV and 4 < y < 5.5 (assuming 7r mass). The events used are 
the 5% most central of all events, with centrality measured using a zero degree calorimeter. 
The products from weak decays such as A's and K 0, s, were only partially rejected with 
approximately 60% rejection efficiency The errors are statistical only. 



how well our resonance gas describes the width of the (almost) Gaussian distribution. From 
the data, we have 

v 2 (N) 

e ^ ; = 2.008 ± 0.009 , (41) 

which we should compare to our resonance gas prediction of 1.5.Q We therefore conclude 
that about 75% of the observed fluctuation is thermodynamic in origin. The contamination 
introduced into the data by fluctuations in centrality could be reduced by analyzing data 
samples with more or less restrictive cuts but the same (N) , and extrapolating to a limit in 
which the cut is extremely restrictive. This could be done using cuts centered at any cen- 
trality. In addition to fluctuations in centrality, there is another experimental (as opposed to 
thermodynamic) factor which could affect the agreement between resonance gas predictions 
and the observed fluctuations. The increase in the fluctuations due to resonances can only 
be detected provided the detector acceptance is large enough to ensure the detection of all 
(or most) of the decay products. NA49 seem to have coverage wide enough to satisfy this 
criterion and a quantitative estimate of losses on its boundaries can easily be made. Our 
resonance gas model predicts that as the centrality cut is tightened, the ratio v 2 be (N)/(N) 
should decrease toward a limit near 1.5. 

Although further work is certainly required, it is already apparent that the bulk of the 
multiplicity fluctuations observed in the data are thermodynamic in origin. Because the 
multiplicity fluctuations are sensitive to impact parameter fluctuations, it may prove dif- 
ficult to explain their magnitude with greater precision even in future. However, the fact 
that they are largely thermodynamic in origin suggests that the effects present near the 

11 In the NA49 data of Table 2, all charged particles are counted whereas we have done our calculations 
assuming that only the charged pions are observed. In our resonance gas model and in the data0, about 
80% of the charged particles in the final state are pions. If we redo the calculation (^), but this time 
define d l r as the number of charged particles (pions, kaons, protons) produced in the i'th decay of the r'th 
resonance, we find that that v^ he (N) / (N) increases, but only by a few percent. NA49 has demonstrated 
that it can study particle identification event-by-event and it may therefore be possible to analyze data on 
charged pion multiplicity in future. 
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critical point, which we study in Sections 5 and 6, could result in a significant nonmonotonic 
enhancement of the multiplicity fluctuations. This would be of interest whether or not the 
noncritical fluctuations on top of which the nonmonotonic variation occurs are understood 
with precision. 

Now, we proceed to pr fluctuations. As we explain in the Appendix, in order to be sure 
that F — 1 when there are no correlations between pions, care must be taken in constructing 
an estimator for w e be(PT) using a finite sample of events, each of which has finite multiplicity. 
The appropriate prescription ( |145|) is to weight events in the event-by-event average by their 
multiplicity. This has not been done in Table 2. However, we show in ( p.47| ) that we can use 
(N) and v e b e (N) to change v e be(pr) as required, and the result is 



to our prediction ( f4"0"|) of 0.68. 

We see that the major part of the observed fluctuation in p? is accounted for by the 
thermodynamic fluctuations we have considered. A large part of the discrepancy is in 
our prediction for the variance of the inclusive single-particle distribution Vi nc (px)- Our 
Vincipr) / {pt) = 0.67 is about 10% lower than that in the data.Q First, this suggests that 
there may be a small nonthermodynamic contribution to the pr-fiuctuations, for example 
from fluctuations in the impact parameter The other source of this discrepancy is the 
blue shift approximation. We have applied a blue shift factor such that (pt) increases from 
281 MeV in Table 1 to 377 MeV as in the data, and in so doing have obtained a value for 
v- inc {pT) which is low by 10%. This may be a reasonable estimate for the error which we have 
introduced by using the blue shift approximation rather than a more sophisticated treatment 
of the effects of flow on the spectrum, which we leave to future work. Such a treatment is 
necessary before we can estimate how much of the 10% discrepancy is introduced by the blue 
shift approximation. Future work on the experimental side (varying the centrality cut) could 
lead to an estimate of how much of the discrepancy is due to impact parameter fluctuations. 

We have gone as far as we will go in this paper in our quest to understand the ther- 
modynamic origins of the width of the inclusive single particle distribution. Another very 

12 As we have already noted, all charged particles are included in the data whereas we have calculated the 
fluctuations for the charged pions alone. We have checked that including the protons and charged kaons 
from the resonance gas increases our prediction for v- mc (pT) / {pt) in the rest frame by only a few percent. 
This small increase in the ratio is likely further reduced once the flow-induced increase in (pt) for the kaons 
and protons is taken into account. Although it would be good to remove this uncertainty completely by 
analyzing a data sample of pions alone, it is already clear that this is not the explanation for the present 
10% discrepancy. 

13 We expect that the fluctuations of an intensive quantity like pt are less sensitive to impact parameter 
fluctuations than are those of the multiplicity, and this seems to be borne out by the data. 



^ebe(PT) = (17.27 ±0.03MeV) 




(42) 



We use this henceforth. We must now compare 



(jV) 1/2 ^ebe(pr) 

(Pr> 



0.751 ±0.001 



(43) 
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important feature in the data is the value of the ratio of the scaled event-by-event variation 
to the variance of the inclusive distribution: 

VF = (A ° 1/Vb \ (PT) = 1.002 ± 0.002. (44) 

The difference between the scaled event-by-event variance and the variance of the inclusive 
distribution is less than a percent .0 This is a remarkable fact, since the contribution of 
the Bose enhancement (see Section |3.1|) to this difference is almost an order of magnitude 
bigger (-y/Fg — 1 is a few percent). Therefore, there must be some mechanism at work which 
compensates for the Bose enhancement. One possible mechanism is the effect of the two- 
track resolution, diminishing the observed number of pions with very similar momenta^]. 
This reduces the ratio F, and NA49 estimates that it is of comparable magnitude to the Bose 
enhancement effect but with the opposite sign. We do not attempt to include either this 
effect or the effect of final state Coulomb interactions between charged pions in our analysis, 
leaving that to the experimentalists. However, we point out that in the next section we 
find another possible origin of this effect. We shall see that anticorrelations due to energy 
conservation and thermal contact between the observed pions and the rest of the system 
reduce F, as long as the system does not freeze out near the critical point. 

In summary, we have shown in this section that the assumption that the system is a 
thermal resonance gas at freeze-out is in reasonable agreement with the magnitude of the 
observed event-by-event fluctuations in the pion multiplicity and mean p?. We will see in 
Section 4 that the effects of energy conservation bring our prediction for y/~F into even better 
agreement with the data. Of course, a number of issues we have touched upon need further 
study: it cannot be otherwise for the first quantitative study of a new set of phenomena. 
The situation is, however, very encouraging. First, RHIC detectors are very well-suited to 
measurements of the fluctuations we have analyzed. Second, although some of the interesting 
effects are small (Fb, for example) with millions of recorded events all Gaussian widths can be 
measured to much better statistical accuracy than even the smallest of the systematic effects 
we have discussed, and will discuss in later sections. Third, the interesting systematic effects 
can be studied by varying the cuts made on the data. For example, considering only low- 
momentum pions one should find the effect of both the resonances and the Bose enhancement 
(see Fig. (]!])) to be several times higher. Also, p n and therefore the Bose enhancement factor 
Fb may be somewhat larger in central events. Detailed study (varying centrality; varying 
cuts in px) may allow experimentalists to separate the effects of Bose enhancement from 
other effects we have described, and will describe later in this paper. Fourth, one can 
significantly widen the types of fluctuations which are analyzed. For example, one can study 
new correlators like the event-by-event cross correlation between multiplicity and mean p?, 
(ApxAN) . We saw in Section 2 that such cross correlation results only from nontrivial 
effects. Finally, it is important to note that we do not expect any of the effects we have 
analyzed in this Section to change significantly near the critical point. 

14 As noted above, because Vi nc (pT) is scaled by the blue shift introduced by the expansion velocity, so 
is <J> PT . This makes $ PT harder to predict than F. However, for convenience, we note that if one uses the 
experimental value of Vi nc (pT), a value y/F = 1.01 corresponds to $ Pr = 2.82 MeV, and the y/~F in the data 
© corresponds to $ Pr = 0.6 ± 0.6 MeV. 
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Our analysis demonstrates that the observed fluctuations are broadly consistent with 
thermodynamic expectations, and therefore raises the possibility of large effects when con- 
trol parameters are changed in such a way that thermodynamic properties are changed 
significantly, as at a critical point. The smallness of the statistical errors in the data also 
highlights the possibility that many of the interesting systematic effects we analyze in this 
paper will be accessible to detailed study as control parameters are varied. 

4 Using an Ideal Gas of Pions as a Thermometer 

To this point, we have assumed that the system does not freeze out close to the critical 
point, and can be approximated at freeze-out as a non-interacting ideal resonance gas. In 
this Section, we take a first step towards understanding how the physics characteristic of 
the vicinity of the critical point affects the event-by-event fluctuations. Along the way, we 
quantify the effects of energy conservation on the pr-fluctuations. This leads to a small 
reduction in \fF far from the critical point, which may be required by the data (JH|). In this 
Section, we consider only the "direct pions" , and as before we treat them as an ideal Bose gas 
at freeze-out. We further imagine that the pions are in thermal contact with the "rest of the 
system" , which is not directly observed and which need not be ideal. The rest of the system 
includes the neutral pions, the resonances, the pions not in the experimental acceptance and, 
most important, the order parameter or sigma field. If freeze-out occurs in the vicinity of the 
critical point, the thermodynamic properties of the sigma field (and therefore of "the rest of 
the system" ) are singular. In the analysis of this section, we imagine that the observed pions 
are an ideal gas even for freeze-out in the vicinity of the critical point, while the equation 
of state and susceptibilities of the rest of the system become singular there. Some of the 
universal critical indices characterizing this singularity are discussed in |]J . The question we 
ask here is how the fluctuations of the pions are affected by being in thermal contact with 
the rest of the system, particularly when the susceptibilities characterizing the sigma field 
diverge. 

A reader who is used to thinking about the 0(4) second order transition may be concerned 
that we are treating the pions and the sigma field so differently. The point is that near the 
critical endpoint which we wish to analyze (and which may occur in nature) the pions 
and sigma are different. The pions remain massive, while the sigma mass vanishes and 
the long wavelength modes of the sigma field undergo critical fluctuations and are almost 
classical. The divergence of the specific heat of the system as a whole is primarily due to 
the fluctuations of the sigma field. The analysis of this section is therefore a reasonable first 
step. What it leaves out, of course, is the fact that the pions, although not massless, do 
interact strongly with the sigma field and are therefore not an ideal gas. We are neglecting 
the direct effects of the pion-sigma coupling. Once these are included, it is not possible to 
make a clean separation between ideal pions and singular rest of the system. We analyze 
the consequences of the pion-sigma coupling in Section 5. 
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4.1 Thermometers, Temperature Fluctuations and Heat Capacity 

Let us step back and recall the text-book formalism describing the measurement of temper- 
ature. A thermometer should be a simple system which has been already calibrated, in the 
sense that we can relate its total energy E to its temperature T via a function T(E) which 
we already know. Instead of E, we could also use any other mechanical observable, like for 
example the volume of the liquid in a liquid thermometer. An ideal pion gas makes a very 
good thermometer because it is a simple system with a known equation of state. Having the 
equation of state, we read off T by measuring a mechanical observable, such as E. 

If the mechanical observable fluctuates, so will the measured temperature. In particular, if 
we measure the total energy, which fluctuates as ((AE) 2 ) = T 2 Cy in the canonical ensemble 
the temperature T(E) will also fluctuate, with ((AT) 2 ) given by (P. Note that given that 
we measure a mechanical observable E rather than T, in order to find the result ([I]) we 
must know T(E). This is possible for a thermometer like an ideal pion gas, but may not be 
possible for the system one wishes to study using the thermometer. One of the questions we 
address in this section is when we measure the energy of the thermometer only (instead of 
measuring the energy of the whole system) which Cy is relevant: that of the thermometer, 
that of the rest of the system, or a combination. 

Suppose now we use our thermometer B to measure the temperature of another system 
A. The measurement consists of bringing the two systems in thermal contact. If the resulting 
system A + B is closed, thermal equilibrium will result. By ergodicity, the thermodynamic 
ensemble will consist of all the states with the same energy, Ea+b, taken with equal proba- 
bility weight. Although the total energy does not fluctuate, the energies of the subsystems 
E A and E B do, subject to a constraint Ea + E b = Ea+b- The probability that the subsystem 
B has energy E B is proportional to the number of states, T, of the system B with energy 
E B times the number of states of A with energy Ea = Ea+b — E B : 

T a+ b{Ea+b) = E r A (E A+ B - Eb)Tb(Eb). (45) 

E B 

Both Y a,b are exponentially growing functions of their arguments (and also the size of the 
system) and their product on the right-hand side of Q4"5|) has a sharp maximum at some value 
of E B . Introducing the entropy S as S(E) = InT(E), we can write for the value of E B at 
the maximum: 

= JL ISa(Ea+b - E B ) + S Bi E B)] = -i§± + H = -i- + (46) 

since the temperature is, by definition, 1/T = dS/dE. We recover the text-book result that 
the temperatures of two systems in equilibrium are equal. Measuring E B , and using the 
known function T B (E B ), we find this (common) temperature. 

Of course, it is not necessary that the system A + B is rigorously closed. In practice it is 
sufficient that the rate of the thermal equilibration between A and B is faster than the rate 
of thermal equilibration of A + B with the environment. 

So far we have only discussed the mean value of E B and, consequently, the mean tem- 
perature. The size of the fluctuations of E B is given by the width of the maximum in 
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Ta(E a+b — E B )Y B (E B ). We need the second derivative: 



M [s ^° ~ Eb) + Sb{Eb)] = Sf + § = i-k + i) ■ (47) 

where C a ,b are the heat capacities of the systems A and B. Thus we find for AE = AE B = 
-AE A : 

((AE) 2 )=T 2 (^- + i-)^ . (48) 

The importance of the result ( f4~8| ) is that the thermometer B allows us not only to measure 
the temperature of the system A, but also the heat capacity of the system A^ In order to 
make such a measurement, we must watch the fluctuations of E B in addition to (E B ). 

Another consequence of (|48D is that when Ca >Cfiwe recover the result for the canonical 
ensemble (|I|). What is important is that the heat capacity Cy appearing in ([I]) in this case 
is that of the thermometer itself, C B , and not that of the measured system, Ca- 

Now, suppose that the system A has a thermodynamic singularity at some temperature, 
as a result of which Ca — * oo. This is precisely the situation which arises near the critical 
point in the idealization of this section: the ideal pion thermometer B is in thermal contact 
with a system A with divergent susceptibilities. Equation (ff8 ) tells us that the fluctuations 



of the energy, which are equal in A and B due to the conservation of energy in A + B, 
will increase as we approach the critical point where Ca diverges. What happens to the 
temperature fluctuations? Remember, that we do not measure the temperature directly, but 
use the equation of state T(E) to read it off from the value of E. If we used the equation 
of state of the system A, T A (E A ), the fluctuations of Ta would decrease and vanish at the 
critical point as discussed in JTJ] , because Ca = dEA/dTA = oo. However, the equation of 
state of the system A is not known to us. Indeed, we are trying to learn about it doing our 
measurements. If we instead use the equation of state of the thermometer T B (E B ) which is 
nonsingular, the fluctuations of T B will increase because the fluctuations of E do, and will 
approach the value determined by (0) with Cy = C B . 

Note that the temperatures of the systems A and B determined through their respective 
equations of state are different on the event-by-event basis. This is not in contradiction with 
thermodynamics which only requires the mean values to agree as in fl-RpD . 

Returning to our idealized system at freeze-out we want to use pions we observe as a 
thermometer, B. The rest of the system, which includes all the other particles, including 
pions not ending up in our detectors, we consider as system A. The singularity of the 
heat capacity occurs in C A , while C B is the heat capacity of the ideal gas and is regular. 
Nevertheless this singularity affects the fluctuations of the pions through (fHf ). The effect of 
the singularity in Ca is an increase in the fluctuation ((AE) 2 ). If one were able to use Ta(E) 
to define T, one would find that fluctuations in T A would decrease at the critical point. 
Using T B (E), or any practical definition of a temperature, leads to fluctuations in T which, 
like those in E, increase. Since what we measure is always a mechanical thermodynamic 
observable, like the total energy E, or the energy per particle, or the transverse momentum 
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One example of such a measurement in a simple lattice system can be found in |27 



25 



per particle, etc., it is not in fact necessary to do a translation to the temperature variable to 
detect a singularity. It is easier to look directly at the fluctuations of observable quantities. 
To this we now turn. 



4.2 The Microscopic Correlator 



As discussed in Section |2.1| , the mean square variations of thermodynamic observables in 
the pion gas are determined by the microscopic correlator: (An p Ank). Once we find this 
correlator we can then use it to calculate any fluctuations of interest. 

For the case of the canonical ensemble this correlator is given by (0) which leads to: 

<(A£) 2 > = £e> p 2 = T 2 C B , ( 49 ) 
v 

as in ([ToD and ([HI) . This corresponds to the case Ca = oo, where B is a (grand) canonical 
ensemble. In the case when Ca is finite (pf) tells us that the correlator (An p Arik) should 
change. A simple derivation of this correlator given below yields the result 

< w = ~ i*c1+%4 ■ m 

This result is easy to understand intuitively and it passes many nontrivial checks. When 
Ca S> Cb = J2 P e p v l/T 2 the second term in ( |50D is negligible and we recover (0). On the 
other hand, when Ca = 0, the system B is closed and the total energy E = J2 P e P n p cannot 
fluctuate. Accordingly, J2 P £ P {An p Arik) = in this case. Note that the correlation is negative 
as it should be, since finiteness of Ca suppresses fluctuations of E, which means that if one 
rip increases, others are more likely to decrease. This negative correlation is therefore a direct 
consequence of energy conservation, and should persist even in systems which are less ideal 
than the one we are analyzing in this Section. 

The microscopic correlator floPP determines the fluctuations of many observables. For 
example, by convolving it with e p 6k as in (|9]) one can derive the result fl48|). This is yet 
another check of fl5"0f ). Note also that the correlation term in (|50|) is down by a factor of 1/V 
(since Ca,b ~ V), where V is the volume of the system. This is also easy to understand: the 
restriction on some linear combination of rip's imposed by energy conservation affects each 
individual n p little if the number of n p 's (i.e., the size of the system) is large. However, the 
contribution of this term to fluctuations of extensive or cumulative quantities is not small, as 
equation (fj8|) shows. This is due to the absence of the Kronecker delta in the second term. 

We now turn to the derivation of the result (|50D . The uncorrelated fluctuations given by 
formula ([/]) follow from the factorizable probability distribution^ 



dP(n p ) = J[dn p exp 1 -— (An p ) 
v I lv v 



(51) 



16 A careful reader may note that eq. literally implies that An p ~ V°, as far as the thermodynamic 
limit 1/V power counting is concerned. n p is also of order V°. If it were the case that An p ~ n p /\/V, 
our assumption that the fluctuations of the occupation numbers n p are Gaussian would be immediately 
justified. Instead, the fluctuations of the occupation numbers n p are not necessarily Gaussian. This can 
be cured by considering, in place of n p , the sum of occupation numbers of a set of modes in a cell (Ap) 3 
centered at p in momentum space, where Ap is fixed as V — > oo. Since the number of modes in such a set 
is (Ap) 3 V — 0{V) and the modes fluctuate independently the central limit theorem will apply and make 
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The energy E = J2 P e p n p in such an ensemble fluctuates according to ((AE) 2 ) = J2 P v l e l — 
T 2 C B . 

Now, if we bring this system into thermal contact with the system A, according to 
( |I^ , |46| , |47|) the probability receives an additional factor: exp[— (AE) 2 / (2T 2 Ca)]- For example, 
if Ca = it becomes a delta-function, meaning that the system B is closed itself, and the 
energy cannot fluctuate. So, we write: 

dP{n p ) = (j[ dn)j exp j - ^ ^ ( An p f j exp j - ^ e p An p ^ 

= (n^ P )/ rfAex p{-E^( A ^) 2 + A E e P A ^} ex p{ T2 ^y^ (52) 

where we have introduced a Lagrange multiplier A. The integration over A should be done 
along the imaginary axis for convergence. 
Completing the squares we find: 



dP(n p , A) = dX drij}j exp j - 7^ ( An v ~ ^ v I € p 

xexpjjW^ + E^y}- (53) 

Now we see that: 

[An p - A^e p ) (An fe - Xv 2 k t^ = v 2 p 5 pk , (54) 



M — ^C^W] 1 - * (55) 



and 

(x(An k -Xv 2 k e k )) = , (56) 
from which we find: 

(An p An k ) = v 2 5 pk + (XAn p )v 2 k e k + (XAn k )v 2 e p - {X 2 )v 2 p e p v 2 k e k 

2f , / \ 2\ 2 2 

= v p b pk + (A )v p e p v k e k 



v p 5 pk - — — — — . (57) 



Q.E.D. 

Now, armed with the equation (|57J|), we can calculate all other fluctuations in our ideal 
Bose gas B in contact with the system A. 



fluctuations of such "smeared" n p Gaussian. Practically, we always convolve n p with a smooth function of p. 
Instead of displaying the smearing of n p explicitly in our notation, we can instead just treat the fluctuations 
of n p as if they are Gaussian, because this will not affect any of the quantities calculated by convolving n p 
with a smooth function of p. 
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4.3 Application: Fluctuations of Mean p T 

As an example of the application of the formula for the microscopic correlator ([50]) we analyze 
the fluctuations of an intensive variable in the ideal Bose gas of pions which we denote q, 
q = Q/N. We wish to see how the fluctuations of q are influenced by the fact that the pions 
are in thermal contact with a system with (possibly singular) heat capacity C A - We shall be 
interested in a particular case where q is the mean transverse momentum p?, but shall use 
the more general notation both for later convenience and to make contact with Section 2 in 
which we discussed e = E/N, another possible q. 

Starting from an equation similar to (|17|), averaging, and using the correlator (|50"D instead 
of (0) we obtain: 

<(Ag) 2 > = 




£^(<z P -<<z» ? ' 



(iV) 2 V p ^ p ^" T\C A + C B ) 



Y. v l e p ( y%- (<?>) 



2 



(51 



The first term on the right-hand side is the same as in (^) with q = e p . This is the 
main contribution to ((Aq) 2 ). We have seen that these thermodynamic fluctuations can 
be described using the variance of the inclusive single-particle distribution and the Bose 
enhancement factor. The second, negative, term in flSBp is the effect of the anti-correlation 
(second term in (|50|) ) induced by energy conservation and thermal contact with the system 
A. This term would be nonzero even if Ca were zero. In this case, it would describe the 
effects of energy conservation on the fluctuations of q in the system B. Thermal contact 
with the system A reduces this term, but it remains important as long as Ca is comparable 
to Cb- It vanishes at the point where Ca diverges. 

In a heavy ion collision, the heat capacity of the pion gas Cb is a sizable fraction of the 
total heat capacity Ca + Cb- The effect fl5"8|) can therefore decrease the fluctuations of px, 
countering the Bose enhancement. This effect will be reduced as we approach the critical 
point where Ca diverges. This will lead to an increase in the event-by-event fluctuations of 
Pt as compared to the variance of the inclusive single-particle spectrum. 

To make the comparison with the Bose enhancement effect estimated in Section |3J] easier, 
we shall express the strength of the effect of the thermal contact in terms of the ratio, Ft 
of the whole expression in curly brackets in fl58|) to the first term in this expression. For the 
fluctuation of mean transverse momentum per event, i.e., for q = p? we find 

0.12 , , 

. (59) 



C A /C 



B I" -L 

for T = 120 MeV and fi n = 0. We see that the effects of energy conservation and thermal 
contact on the fluctuations of an intensive quantity like pt are smaller than the effects 
on the fluctuations of the energy in (f48|). Several obstacles make it difficult to use ( |59"D 
quantitatively. First, some dilution of the effect is to be expected because less than half of 
the pions which are observed are direct. Second, it is a little bit difficult to know how to 
estimate Ca/Cb, because we have analyzed such an idealized situation. The system A should 
certainly include the neutral pions in the same region of rapidity as the observed charged 
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pions; however, should the pions at different rapidities be included? The total number of 
pions, neutral and charged, in a central event at the SPS at 160 AGeV is about ten times 
larger than the number of charged pions per event used in NA49's present analysis. This 
suggests that Ca/Cb is at the very most 10. If we take Ca/Cb ~ 3 for orientation, which 
can be justified if one assumes that Ca includes the heat capacity of the resonances and that 
of the neutral pions in the same region of phase space as the observed pions, we find Ft — 1 
of the order of —3%, before taking into account the dilution by non-direct pions. The effect 
is comparable in magnitude to the Bose enhancement, acts in the opposite direction, and 
should be reduced near the critical point at which Ca diverges. 

A divergent specific heat is only possible in an infinite system. In Section 5 we will 
estimate that in a realistic heavy ion collision, finite size effects suggest that near the critical 
point the sigma contribution to Ca could be as much as a factor of ~ 6 2 larger than the 
contribution of a typical light degree of freedom. This suggests that Ca could easily increase 
by as much as an order of magnitude at the critical point, reducing the anti-correlation in 
(AripArik) and the negative contribution to Ft by the same factor. 

The effects of thermal contact can be distinguished from other effects, like those of finite 
two-track resolution, which also counter the Bose enhancement effect because of the specific 
form of the microscopic correlator floTf). The effect of energy conservation and thermal 
contact introduces an off- diagonal (in p k space) anti-correlation. Although our estimate 
of the magnitude of the effect suffers from a variety of uncertainties introduced by the 
idealizations used throughout this Section, the existence of this off-diagonal anti-correlation 
is robust. It arises simply due to energy conservation: when one n v fluctuates up others 
must fluctuate downward, and it is therefore more likely that nk fluctuates downward. If Ca 
is increased, then the system A can more easily supply the energy needed for the upward 
fluctuation in n p , and the anti-correlation between n p and nk is reduced. Preliminary analysis 
by NA49 suggests that some amount of such anti-correlation is observed in the data[^£|. It 



will be interesting to compare the magnitude of any effect observed in the data with our 
estimates. We leave this to future work. If it is possible to separate this effect from other 
effects because it is an off-diagonal anti-correlation, then a measurement of this effect would 
yield an estimate for the effective value of the ratio Ca/Cb at freeze-out. 

Note that Ft increases near the critical point, but it increases towards a finite value 
(namely 1.) In contrast, in Section 5 we will explore effects which result in the divergence of 
an analogously defined ratio F a at the critical point. 



4.4 Another Application: T 

In this subsection we introduce another measure of the temperature of the pion gas, T. Our 
new variable T is well-defined on a single event, and has the property that (T) is related 
to the slope parameter. We have found that although (pt) is related to T, the fluctuations 
((Apt) 2 ) are not at all like the fluctuations ((AT) 2 ) in (0). We now show that ((AT) 2 ) 
also does not behave quite like ((AT) 2 ). The reader should expect this, since we argued on 
general grounds that (|l|) can only be obtained from a mechanical observable if the equation 
of state of the system A is known. Still, it is nice to confirm this using an example of an 
observable which is a less straightforward intensive quantity than just e = E/N or pt- 
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We define for each member of the ensemble (i.e., for each event) independently: 



(60) 



where 



"SCO 



oe p /T _ I ■ 



(61) 



and dp is some function of p which we can choose for convenience later. Then for each event 
individually we can define a temperature, T, which is found by minimizing x 2 {T) for this 
event: 



dx\T) 
dT 



0. 



(62) 



T=T 



It is clear that mean value of T over all events for the ideal Bose gas will coincide with the 
actual temperature of the gas T. But, since T is defined for a single event, rather then for 
the whole ensemble, it fluctuates! 

As before, the fluctuations in T are determined by fluctuations of n p . For small fluctuation 
AT we can write: 







A 



dx\T) 
dT 



A 



T=T 



/ 9n° n 

Y AT— ^ 

V V 9T 



- An r 



dn° 1 

£K(T)-V p 

\ dn° p 1 



8Toj 



dT o* 



(63) 



where we omit the tilde on T. Note that even if contains dependence on T, it will be 
always multiplied by (n p — n°) which is zero to the relevant order in the size of the fluctuation. 
Now we need: 



dn° p (T) 
dT 



—e p n° p (l + n° p ) 



rp2 e P V V 



according to (|). We rewrite ( |B"3"D as: 



AT 

J^2 



(64) 



(65) 



We can carry on with an arbitrary a p , but let us make the following choice: a p = v p . 
This choice makes a lot of sense if one recalls that in the standard definition of \ 2 one divides 
each square deviation term by its normal square deviation (which is usually obtained from 
experimental error, and which here we know to be v p from the fluctuations of n p ). This 
choice simplifies formulas. 



Now we square (65[), average over events, and restore the tilde on T: 

2 

= £ P £k (An p An k ) , 



((AT) 2 ) 



E^ 2 



(66) 
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which, according to (|50|) means: 
((AT) 2 ) 1 C A 



T 2 C B C A + C_ 



B 



(67) 



We see that fluctuations of T, like those of pt, increase towards the critical point of the 
system A, where C A — * oo, approaching a finite constant. When C A is infinite, the system 
B (the Bose gas) is in the canonical ensemble, and the fluctuations of T are given precisely 
by (0), with the specific heat Cb in the denominator.^ 

Different definitions of the temperature T can be devised (using different choices of ex 2 ). 
They will lead to different temperatures for a given event, which are the same in the mean 
(and equal to T), but different in the size of their fluctuations. All these fluctuations will 
increase somewhat at the critical point, but will not diverge there as they are controlled 
there by Cb (or some other property of the thermometer B) which is nonsingular. 

4.5 Two Further Applications: ((AiV) 2 } and (AN/Ap T ) 

Once we understand how some physical effect influences the microscopic correlator (An p Ank), 
we can calculate the fluctuations of many different observables. The task then is to look for 
observables in which the effect of interest is large, and which are of practically utility in the 
sense that they are easily accessible to experimental analysis. We give two further simple 
examples here. 

In Section 4.3, we analyzed the fluctuations of an intensive quantity, p T , and obtained 
the expression (pffi ) for ((Apt) 2 )- Similarly, we can use the microscopic correlator ( pT|) to 
analyze the fluctuations of the extensive quantity N, and obtain 



<(AA0 2 > = £ 



T*(C A + C B 



12 



L V 



(68) 



The first term is the ideal Bose gas result, and the second term is the correction due to 
thermal contact and energy conservation. For T = 120 MeV, /x w = 0, the effect of the second 
term is to multiply ((AiV) 2 ) by a factor of [1 - 0.20/(1 + C A /C B )]. Note, however, that the 
multiplicity fluctuations of the pions obtain from the resonance gas which we analyzed in 
Section 3 are dominated by the pions from those resonances which decay into more than one 
pion. Doubling the contribution of the direct pions to ((AiV) 2 ) in the calculation ( |3"2"D only 
increases ((AiV) 2 ) by 10%. The effect of thermal contact and energy conservation on the 
direct pions seen in (|68|) is therefore a very small contribution to the total ((AiV) 2 ) of (|32|). 

We saw at the end of Section 2 that cross correlations between intensive observables 
and N are of interest, because they vanish in a classical ideal gas. We therefore use the 
microscopic correlator fl50|) to calculate 



<^> = M{5**-< 8>) -nCT 



p 

■ p 



v l^p{% - (?}) 



.(69) 



17 The heat capacity Cy in this case is that at fixed fi. It is remarkable that if we define t by simultaneously 
fitting two variables in (|60|), T and ft,, the resulting T will again fluctuate according to ([l]), but with heat 
capacity at constant N. We leave this as an instructive exercise for our reader. 
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For q = pt, T = 120 MeV, /i, = 0we find 

< A f ApT> - -0-021 fi + * a r ) ■ (-0) 

(pr) V Ca/Cb + IJ 

As we saw in Section 2, correlations like this arise only due to nontrivial effects, and are 
generally small. In this case, we see that (for Ca ~ 3Cb) the effect of energy conservation 
and thermal contact is ~ 2.5 times as large as that due to Bose enhancement. This suggests 
that this correlation would be a very interesting quantity to use to look for the critical point. 
It is small in magnitude, but even after the dilution of the direct pions by those produced 
in resonance decays are taken into account, it may change by a large factor near the critical 
point where Ca — > 0. 

In conclusion, the effects of thermal contact and energy conservation on the pions could 
either be found directly, by detecting the anti-correlation in the microscopic correlator 
(An p An k ). Or, the resulting effects on ((Ap T ) 2 ), ((AT) 2 ), ((AiV) 2 ) or (Ap T AN) which 
we have estimated may be discovered, likely by seeing them change as control parameters 
are varied. 



5 Pions Near the Critical Point: Interaction with the 
Sigma Field 

In the previous section, we made the assumption that the "direct pions" at freeze-out could be 
described as an ideal Bose gas. We do not expect this to be a good approximation if the freeze- 
out point is near the critical point. The sigma field is the order parameter for the transition 
and near the critical point it therefore develops large critical long wavelength fluctuations. 
These fluctuations are responsible for singularities in thermodynamic quantities. In the 
previous section, we analyzed this situation by pretending that the only effect on the pions 
was due to thermal contact with a heat bath with divergent susceptibilities. In this Section 
we take the next logical step, and consider the effect of the classical critical fluctuations 
on the pions through the aim coupling. It would be strange if, as in the previous section, 
the properties of the pions remained regular in the thermodynamic limit in the presence of 
the nonanalytic behavior of the sigma field. We will see that the fluctuations of both the 
multiplicity and the mean transverse momentum of the pions do in fact diverge at the critical 
point. 

We then estimate the size of the effects in a heavy ion collision. This requires first 
estimating the strength of the coupling constant G, and then taking into account the finite 
size of the system and the finite time during which the long wavelength fluctuations can 
develop. The pion fluctuations induced by the Goixix interaction are divergent and are 
therefore the dominant fluctuations in an infinite system. In the finite system of interest, 
we find that the momentum fluctuations are large enough to be easily detectable, but not so 
large as to seriously jeopardize our treatment, which considers the effects of the interaction 
only to lowest order. It is for this reason that we have first analyzed all effects other than 
those introduced by the Garni interaction, and now add these effects in. The multiplicity 
fluctuations are large enough that in this case a treatment which goes beyond lowest order 
in G seems called for. We leave this to the future. 
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5.1 Microscopic Correlator 

As before, we shall derive the expression for the microscopic "master" correlator, (An p Arik) , 
which can then be used to calculate fluctuations of various observables. We neglect the 
effects considered in the previous Section, as they can be added to the effects of this Section 
at the end. We concentrate on the fluctuations of the sigma field, the fluctuations of the 
pion occupation numbers, and the otttt coupling. The long wavelength fluctuations of the 
sigma field which are responsible for the singular effects of interest are classical. 

The effective potential Q determines the probability distribution of the classical field a 
through 

dP(a)=d<rexp{-^}. (71) 

This equation can actually be thought of as the definition of fi(er). The effective potential is 
extensive, but for convenience we set the volume V = 1 in the calculations to follow, although 
we will restore it explicitly in our results. (Note that throughout previous sections, we had 
set V = 1 implicitly. The momentum sum J2 P should always be read as V J d 3 p/ {2tt) 3 .) Let 
us consider small fluctuations of the field a around the minimum of Q(a). We can then 
expand the effective potential Q(a) around a = 0. The first terms will be 

77? 2 rfT 2 

Q(a) = + Ga:n 2 : +0(a 3 ) , (72) 

where we have temporarily omitted terms independent of a (such as m 2 7r 2 /2).Q The second 
term is the interaction between sigmas and pions. The coupling G has the dimensions of 
mass, and its magnitude near the critical point will be estimated below. The notation :: 
signifies tadpole subtraction: : tt 2 := tt 2 — (tt 2 ), which makes sure that the minimum of a is 
not shifted as we shall see below. (The notation tt 2 itself is itself somewhat symbolic, as it 
represents / d 3 x7r(x)7r(x).) Thus we have: 

dP(a) = daexp j-^ - : tt 2 : J (73) 

Now, the field tt also fluctuates. Let us determine the corresponding (joint) probability 
distribution. In the previous section we used the probability distribution for the occupation 
numbers, and we begin by translating the fluctuations of the field tt into fluctuations of the 
occupation numbers. We write, doing the usual Fourier transform: 



TT 



2 £ K\ 2 - (74) 



V 



We can relate the Fourier components tt p to the occupation numbers n p . It is clear that 
n p ~ \k p \ 2 . The coefficients can be determined, for example, by using 

Z = J Vtt exp | - fJ T dt J dV [| (<V) 2 + X -m\ vr 2 ] j = Z T=0 J[ (l - e"^ /T ) . (75) 



18 Clearly, the fluctuations of a are not small. We shall proceed with the assumption that the higher-order 
terms in yield subleading contributions to the singular effect we seek. We shall return to this point in 



Section 5.4. Also note that we consider only the zero momentum mode of the field a. This can be justified 



in a diagrammatic approach, which can also handle nonzero momentum modes of a. 
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Differentiating In Z with respect to m 2 we find 

(vr 2 ) = E-K), (76) 
p u p 

up to the temperature independent vacuum contribution (equal to ^ p u p /2) from lnZr=o- 
So we have 



7T 



2 



1 



E - V (77) 



p ^p 



Note now that (n 2 ) = J2 P ( n p) I luJ p 7^ 0. So, unless we subtract (ti 2 ) the minimum of a will be 
shifted from the origin (this subtraction will also take care of the vacuum fluctuations). We 
have: 

- 2 -(vr 2 ) = E^- (78) 
p ^p 

Now, putting everything together, we find the joint probability distribution for the sigma 
field and for the pion occupation numbers: 



dP{°, n p ) = da (n dr^ exp j- £ ^(An p ) 2 - ^ £ ^ - |^ 2 } . 



(79) 



This is a very important formula which will allow us to calculate the fluctuations. 

The measure dP(a, n p ) is Gaussian, which is very helpful. Completing the squares, we 
find 



dP(a, n p ) = da (j[ dn^j exp i - E g~2 ( A 



n p + — ^ 



«-p{-(#-^)^} <*> 

Before we make the final and the simplest step, let us make two side notes. 

The equation ( |80| ) shows that the interaction with a shifts mean occupation numbers by 

S(n P ) = . (81) 

As the reader might have guessed already, this must be due to the shift of the mass of the 
pions linear in a, which can be seen from ([7^): 

= (82) 
It is trivial to evaluate the change in (n p ) induced by the change of the mass. Since e p = 

p2 _|_ m 2 _ ^ we f| nc l : 

5(n p ) = -K)(l + (n p )) 6 -^ = (83) 
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which is the same as (|81|) with (|82|). The fluctuations of the sigma field will have further 
affects on the pion occupation numbers, but these are higher order in G and we neglect them 
here. 

The second side note is more important. We see from (|80D that the mass of the a field is 
corrected by the fluctuations of the pions: 

~ 2 _ 2 V p 



m„ = m„ — 



E-4- (84) 



T 2 „ ujI 



Diagrammatically, this corresponds to the thermal one- loop diagram a — > tttt — > a. The 
physical mass of the sigma is rh a , to the order in which we are working. This is the mass 
which vanishes at the critical point. We shall omit the tilde in the following. 

Finally we can read off the following expectation values from the probability distribution 

IUn r + ^i\U nt + ^l\\^ liSrk] (85) 



T u p J \ T iu k/ 



<^ 2 > = ^; (se) 



(aAn p ) = -(a 2 £i. (87) 
Tu p 

This gives 

1 G 2 vM 



(An p An k ) = v p 5 pk + 



ml T uj p uj k 



We see that the coupling of the pions to the sigma field leads to a singular contribution to the 
correlator of the pion fluctuations as we approach the critical point at which m a = 0. The 
first term on the right hand side describes the variance of the inclusive distribution and the 
Bose enhancement effect as we saw in Section 3. The additional terms which we discovered 
in Section 4 could now be added to the right hand side. It is of course the new, divergent 
term on which we shall focus our attention. 

One can represent both terms in this equation diagrammatically as in Fig. 0. The singular 
term is due to the exchange of the sigma in the process of forward pion-pion scattering. This 
results in a characteristic 1/m 2 singularity. A different way of deriving the formula for the 
correlator would be to do a straightforward diagrammatic expansion of (An^An^.). This will 
also allow one to include the effects of the nonzero momentum modes of the sigma field. The 
second term in (|88"D is the most singular one in this correlator as m a — > because it involves 
the exchange of the sigma field with zero momentum. We defer an analysis of less singular 
terms using the diagrammatic approach to future work. 
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Figure 2: Diagrammatic representation of the right hand side of the correlator fl8Sj) . The 
crosses represent the insertions of An p =: tt p tt* :. The solid/dashed lines are the pion/sigma 
field thermal propagators. 



5.2 Application: Fluctuations of Mean pr 

Using (|8"5|), we can determine the fluctuations of any thermodynamic observable which can 
be expressed in terms of the pion occupation numbers. For a generic intensive observable 
q = Q/N we find 



<(Ag) 2 > = 




where we have displayed the factor of V explicitly to show that both terms are of the same 
order in V. (Recall that J2 P ~ y V J d 3 p/ (2tt) 3 .) As before, the intensive observable of primary 
interest will be the mean transverse momentum in an event, q = Pt- It is clear from ([89]) 
that the fluctuations of pt increase near the critical point and diverge at the critical point, 
where m a vanishes. We will give a quantitative estimate of the effect in Section |5.5| . We 
must first estimate the size of the coupling G and the value of m a near the critical point .0 



10. 



(89) 



5.3 The Size of the Coupling G 

The strength of the singular contribution to the pion correlator near the critical point depends 
on the size of the coupling G between the a and n: 

C I = Ga7T i 7i i (90) 

where the isospin index i = 1,2,3 is summed. We first make a phenomenological estimate 
of the magnitude of G in vacuum, and then estimate by how much G is reduced near the 
critical point E. 

The value of this coupling in the vacuum can be inferred independently from two con- 
siderations: (i) from the relationship between the sigma and pion masses and /„.; (ii) more 
directly, from the width of the sigma. We shall use both and compare. 

One way of estimating the vacuum value of G is to use the Gell-Mann - Levy linear 
sigma model f29|, in which the Lagrangian describing the dynamics of the four component 
field 4> a = (0o, t) is given by 

C= Jd A x l^fadrfc - j - v 2 ) 2 + Hfo] , (91) 



19 Similarly, one can also calculate the fluctuations of T denned in Section 4.4 and show that these fluctu- 
ations also increase near the critical point and diverge at the critical point. 
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where the 0(4)-breaking field H is proportional to the quark mass: H = m(%pip) / '/„.. The 
vacuum expectation value (<f> ) is nonzero and should be set equal to f n . The a field is then 
defined by a = <p Q — (</> ). Setting f n = 93 MeV, m n = 135 MeV and m a = 600 MeV fixes 
all three parameters in the potential, and in particular yields A = 20.0. If we rewrite (|9l"|) in 
terms of a, we find a term Xf^an 2 , from which we conclude that 



This value for G seems large at first sight, but such a large coupling is in fact required 
by experiment. In order to see this, we evaluate the width of the sigma due to its tree- level 
decay into two pions, and find 



where we have used m a = 600 MeV as above. The width of the sigma is known experi- 
mentally to be so large that this "particle" is only seen as a broad bump in the s-wave tt-tt 
scattering cross-section. An estimate of 300 MeV for this width is therefore reasonable. We 
conclude that the vacuum aw coupling must be at least as large as G ~ 1900 MeV, since 
the sigma would otherwise be too narrow. 

Our estimate makes it clear that the vacuum value of G would not change much if one 
were to take the chiral limit m — > 0. The situation is different at the critical point. Taking 
the quark mass to zero while following the critical endpoint leads one to the tricritical point 
P in the phase diagram for QCD with two massless quarks. At this point, G vanishes as we 
discuss below. This suggests that at E, the coupling G is less than in vacuum. Our goal in 
the remainder of this subsection is to use what we know about physics near the tricritical 
point P to make an estimate of how much the coupling G is reduced at the critical endpoint 
E (with the quark mass m having its physical value), relative to the vacuum value G ~ 1900 
MeV estimated above. 

We begin by recalling some known results. (For details, see Refs. f|], ||, |lj.) In QCD 
with two massless quarks, a spontaneously broken chiral symmetry is restored at finite 
temperature. This transition is likely second order and belongs in the universality class 
of 0(4) magnets in three dimensions. At zero T, various models suggest that the chiral 
symmetry restoration transition at finite fi is first order. Assuming that this is the case, one 
can easily argue that there must be a tricritical point P in the T\x phase diagram, where 
the transition changes from first order (at higher /i than P) to second order (at lower //), 
and such a tricritical point has been found in a variety of models. [0, ||, |50[ The nature of 
this point can be understood by considering the Landau-Ginzburg effective potential for <p a) 
order parameter of chiral symmetry breaking: 



The coefficients a, b and c > are functions of \i and T. The second order phase transition 
line described by a = at b > becomes first order when b changes sign, and the tricritical 
point P is therefore the point at which a = b = 0. The critical properties of this point can be 
inferred from universality |2|, [J, and the exponents are as in the mean field theory (|9~4"D. We 



G = Xf w ~ 1900 MeV . 



(92) 




(93) 



^(0a) = 7T0a0a + T(0a0a) 2 + ^(4>a4>a) 3 ■ 

2 4 6 



(94) 
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will use this below. Most important in the present context is the fact that because (0) = 
at P, there is no cttttt coupling, and G = there. 

In real QCD with nonzero quark masses, the second order phase transition becomes a 
smooth crossover and the tricritical point P becomes E, the second order critical endpoint of a 
first order phase transition line. Whereas at P there are four massless scalar fields undergoing 
critical long wavelength fluctuations, the a is the only field which becomes massless at the 
point E, and the point E is therefore in the Ising universality class|2|, ^J. The pions remain 
massive at E because of the explicit chiral symmetry breaking introduced by the quark mass 
m. Thus, when we discuss physics near E as a function of /i and T, but at fixed m, we will 
use universal scaling relations with exponents from the three dimensional Ising model. Our 
present purpose, however, is to imagine varying m while changing T and \i in such a way as 
to stay at the critical point E, and ask how large G (and m n ) become once m is increased 
from zero (the tricritical point P at which G = m n = 0) to its physical value. For this task, 
we use exponents describing universal physics near P. Applying tricritical scaling relations 
all the way up to a quark mass which is large enough that m n is not small compared to T c 
may introduce some uncertainty into our estimate. 

In order to determine the trajectory of the critical line of Ising critical points E as a 
function of quark mass m,Q it is sufficient to consider the effective potential only as a 
function of the single component 0o = of the four-component order parameter. When the 
quark mass is nonzero we can add terms containing odd powers of 0: 3 and 5 , in addition 
to just 0. We shall assume that the linear term provides the leading effect, and check this 
assumption for self-consistency a posteriori. So, we have at nonzero m 

n(0 = _ m + ^ + ^ + ^ (95) 

We assume that the units of mass are chosen in such a way that the coefficient of the linear 
term in (^) assumes this simple form. That is, instead of writing it as H — mM 2 , we 
are using units with M = 1. Stable or metastable thermodynamic phases are described by 
minima of Q, at which Q' = 0. At the critical point E, Q' — and in addition both Q" and 
Q'" vanish. This is because three roots of the polynomial Q'(4>) coalesce (two minima of Q 
and one maximum in between). So, we have three conditions: 

n' = -m + a(0) + 6(0) 3 + c(0) 5 = 0; (96) 

Q" = a + 36(0) 2 + 5c(0) 4 = 0; (97) 

Q!" = +66(0) + 2Oc(0) 3 = 0. (98) 

These conditions allow us to express a, b and (0) (the value of at the minimum), as 
functions of m and c. We neglect any change in c; it is the vanishing of m at P which is of 
interest to us. Solving these equations by working up from the last to the first and keeping 
only the exponents (neglecting pre-factors) we find: 

a ~ m 4/5 ; -b ~ m 2/5 ; (0) ~ m 1/5 . (99) 

20 See Ref. |5l) for a derivation of the analogous line of Ising points emerging from the tricritical point in 
the QCD phase diagram at zero /iasa function of m and the strange quark mass m s ■ This tricritical point 
can be related to the one we are discussing by varying m s . 
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The power 1/5 is easy to understand: it is 1/5, where 5 = 5 for the 6 potential. 

At m = 0, the tricritical point P at a = b = has ((f)) = 0; the expressions (|99|) describe 
how the location of the critical point E in the ab plane, and the value of ((f)) at E, change 
as m is increased from zero. From these, we will determine how m^ and G at E vary with 
m, after two asides. First, note that from these universal arguments we learn nothing about 
the location of the tricritical point a = b = in the T/x plane. One can only make rather 
crude estimates of the position of this point, as we have done in Ref. M. Our main purpose 
here and in |IJ is to tell experimentalists how to find P, so that they can find it and tell us 
where it is. Second, we must estimate the size of the 3 and 5 terms we have neglected. 
Assuming that both terms come with coefficients which are at least linear in m (higher odd 
powers of m are possible, but will make the size of these terms even smaller) and using the 
m power counting of (0), we see that while all the terms in (pop are (9(m 6//5 ), the 3 and 
5 terms contribute at most (9(m 8//5 ) and 0(m 2 ) respectively. 

To follow m-x and G, we need the full dependence of Q on the 0o and it fields: 

0(0, 7T) = -m4> + °-{<f + 7T 2 ) + \((f? + 7T 2 ) 2 + Utf + 7T 2 ) 3 . (100) 

I 4 o 

For the pion mass, we need to expand around = ((f)) and 7r = and collect order 7r 2 terms: 

fi=^(a + 6(0) 2 + c(0) 4 ) + ... . (101) 

We can now read off the pion mass: 

ml = a + 6(0) 2 + c(0) 4 , (102) 
which, according to (|96|), means: m 2 = mj ((f)). Using ( p9|) we find: 



ml ~ m 4/5 . (103) 

Assuming that the dimensionful factor in this formula is of the same order of magnitude 
as the one in the zero T and \x formula m 2 ~ m we conclude that the pion mass does not 
change much from its vacuum value, and is likely to be very slightly bigger (by a factor of 
order (Aqcd/ 1 ^) 1 ^ 10 )- This is similar to what is known to happen near T c for /x = 0. (See 
Ref. |31| for a review.) 

To determine the constant G, we need to collect the <T7r 2 terms in ( 100|) where, as before, 
a = (f> — ((f)) . Only the last two terms contribute and we find 

n = o-Tr 2 (b((f)) + 2c((f)) 3 ) + ... . (104) 

This means G = b(<p) + 2c(0) 3 which, according to (]98f) , gives G = 2b ((f))/ 5 = — 4c(0) 3 /3. 
Using the m power counting (P5| ) we find 

G ~ m 3/5 . (105) 

Thus the coupling G is suppressed compared to its "natural" vacuum value G vac by a factor 
of order (m/A Q c D ) 3/5 . Taking A QCD ~ 200 MeV, m ~ 10 MeV we obtain our estimate 

G E ~ ~ 300 MeV . (106) 
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The main source of uncertainty in this estimate is our inability to compute the various 
nonuniversal masses which enter the estimate as prefactors in front of the m dependence 
which we have followed. In other words, we do not know the correct value to use for Aqcd 
in the suppression factor which we write as (?ti/Aqcd) 3//5 - 



5.4 Finite Size and Finite Time Effects 

The final ingredient needed for the estimate of the size of the effect described in Sections 5.1 
and 5.2 is an estimate of m a . We found that (An p Arik) is infinite when m a = 0, at E. This 
singularity occurs because the correlation length £ of the sigma field is infinite. In practice, 
however, there are important restrictions on how large £ can become. The fireball created 
in a heavy ion collision has a finite size and lives for a finite time; both restrict £. Similar 
considerations affect the estimate of the size of the effect described in Section [4.2| . There, 
we found an anti-correlation in (An p Arik) which vanishes as the specific heat of the system 
diverges. The limit on £ introduced by finite size and finite time effects also limits how large 
the heat capacity Cy becomes. 

We discuss finite size scaling first. If the system is infinite, a singular thermodynamic 
quantity such as Cy diverges at the critical point. If the system is large relative to microscopic 
scales (~ 1 fm in our case) but finite, then Cy exhibits a peak at the critical point which 
becomes narrower and higher as larger and larger systems are considered. Finite size scaling 
analysis tells us how the magnitude of the peak scales with the system size. The scaling 
postulate tells us that the singular parts of all observables are due to the diverging correlation 
length £ and can be characterized by an appropriate critical index: Q S i ng ~ £ A<3 , where Q 
could be Cy or some other quantity which diverges at E. In a finite system the growth of the 
correlation length is limited by the size of the system: £ max ~ R. Therefore, the magnitude 
of the singularity in a given thermodynamic quantity (the height of the peak) depends on 
the size of the system as 

Q max ~R AQ . (107) 

Similarly, if the system is not allowed enough time to equilibrate, the singularity is also 
smeared. The magnitude of the singularity in this case can be estimated using finite time, or 



dynamic, scaling [32]. In this case the scaling postulate tells us that the typical equilibration 
time diverges at the critical point (critical slowing down), with this divergence related to 
that of the correlation length by t ~ £ A *. Reversing this relation, we conclude that if the 
typical time allowed for the system to equilibrate is limited to r, the correlation length can 
only grow up to £ max ~ r 1 ^'. Thus, in this case 

Qmax ~ r A «/ At . (108) 

The calculation of the numerical prefactors in ( |107| ) and ( |108| ) requires precise knowledge 
of the QCD dynamics and is not feasible at this time. The exponents, however, are universal 
and can be understood by relating them to suitable exponents in the three-dimensional Ising 
model. For example, the exponent for the specific heat Cy at the end point E was determined 
in [|: 

A Cv = (1) « 2. (109) 

WJ 3d-Ising 
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A Cv is not given by the (smaller) exponent (a/z/) 3d „ Ising because of the obliqueness of the 
first order phase transition line relative to the T axis on the phase diagram as explained in [0. 
The idea is that at the critical point, Cy = d 2 Q/dT 2 is related to some linear combination 
of the Ising model susceptibilities d 2 Q/dt 2 , d 2 Q/dtdh, and d 2 Q/dh 2 where the Ising model 
temperature axis t and magnetic field axis h are oblique relative to the T and n axes. Cy is 
controlled by the most divergent of the three Ising model susceptibilities, which is d 2 Q/dh 2 , 
and ( |109p results. 

The dynamic scaling exponent A t , which is often called z, is also universal. The dynamic 
universality class of a system is sensitive to details of the dynamics such as whether the order 
parameter is or is not conserved and whether the system has other conserved quantities. 



The determination of A t is a rather involved problem in some cases [32]. If we assume that 



QCD at the critical point E falls into the dynamic universality class of the gas-liquid phase 
transition (model H in the classification of Hohenberg and Halperin[]32]]) the exponent A t 
can be estimated as A t w 3. It may therefore turn out that because A t > 1 the finite time 
scaling restriction ( |108| ) may be somewhat more restrictive in a heavy ion collision than the 
finite size scaling restriction ( |107| ). 

Let us estimate some typical numbers for central PbPb collisions at the SPS. We start 
with an estimate for the relevant size in the longitudinal and transverse directions beyond 
which £ cannot grow. The longitudinal expansion extends the longitudinal size of the fireball 
considerably, but regions with different rapidities freeze out at different times, and a homo- 
geneous freeze-out at a single freeze-out time for all rapidities is not a good approximation. 
A similar (although not identical) problem has already been faced in two pion interferom- 
etry, which provides sizes (and durations) not of the whole system, but of a "patch" large 
enough that particles emitted from it can still interfere. The size of such patches depends 
somewhat on the direction and magnitude of the total momentum of the pion pair used in 
the interfero metric measurement, but is an approximate measure of the size of the system 
over which freeze-out is homogeneous. The longitudinal size of such a patch for central PbPb 
collisions at 160 AGeV is estimated to be|25] 2R^ ^12 — 14 fm. At its ends, the rapidity 
difference is already about 1. We therefore estimate that sigma field correlation length is 
limited by finite size effects to be less than 2Rl. 

The size Dt in the transverse direction beyond which £ cannot grow can be estimated 
in two ways. The initial size is that of the diameter of the nuclei, -Dpb ~ 14 fm. The 
transverse (radial) flow makes the physical size of the freeze-out surface larger than the 
nuclear radius, by 30-50% at freeze-out. Therefore, it must be the case that D T < 20 fm. 
This is, however, an overestimate. Because of the relativistic transverse expansion, regions 
with different positions in the transverse directions cross the transition region and then 
freeze-out at different times. Therefore, as for Rl above we can use the size of "patches" 
observed via two particle interferometry as a guide, the sigma field correlation length in the 
transverse directions to be less than = 2R? ~ 10 — 12 fm[^]. It therefore seems that 
the relevant longitudinal and transverse length scales at freeze-out are about the same, and 
we conclude that based on finite size restrictions alone 

y < 12 (110) 

where I is the "microscopic length" of order 1 fm. 
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We now turn to the restriction on the correlation length which arises from the fact 
that the matter created in a heavy ion system does not enjoy an infinite period of time in 
which to equilibrate. The expansion time can be defined through the corresponding "Hubble 
constant" 

1 ds , 
H= — = —. Ill 

th sat 

where s is entropy density. We use the entropy density in the definition because the total 
entropy is conserved during adiabatic expansion, and we are assuming that by the time the 
system is traversing the transition region and then freezing out, the expansion can be treated 
as adiabatic. Hydrodynamic models|20| suggest that at SPS energies, heavy ion collisions 
have th ~ 10 — 20 fm/c. If we simply use this value of th and neglect dimensionless factors 
in the scaling relation ( |108| ) we would find 

7<(t) 1/A '~ 2 - 5 - < 112 > 

In spite of the long expansion time, the relatively large value of the dynamical exponent 
A+ can make the finite time restriction more severe than the finite size one ( ITTOD . In other 
words, although the size of the system allows the correlation length to become as large as 
12 fm, there may not be enough time for such long wavelength fluctuations to develop due 
to critical slowing down. The estimate ( |112| ) is suspect for several reasons. First, there 
may be a large dimensionless proportionality constant in this relation which is unknown to 
us. (In contrast, the finite size estimate ( P-1CQ is a consequence of geometry, and unless the 
homogeneous region is larger than we estimate, it is unlikely that the finite size bound on £ 
can be evaded.) Second, in making the estimate ( |112[ ) we have estimated the "bad" effect of 
critical slowing down, namely the fact that a 12 fm correlation length will take longer than 
12 fm to develop, but we have not taken into account a compensating "good" effect of critical 
slowing down: because of the large specific heat, the system will spend an unusually long 
time with a temperature in the vicinity of the critical point. Because of the uncertainties in 
( |112| ), we shall use £ max ~ 6 fm as a rough estimate of the largest correlation length possible 
if control parameters are chosen in such a way that the system freezes out close to the critical 
point. More detailed study of the time evolution of the temperature of the system, of the 
appropriate choice for th , and of the dimensionless factors in ( |112j ) are required in order to 



properly estimate whether finite time effects restrict £ max further. 

Since the thermal contact effect discussed in Section 4 depends on the divergent heat 
capacity Ca, we need to estimate how large Ca can get, given the finiteness of the system. 
Using the exponent ( |109| ), we can estimate the ratio of the maximum value of that part of 
Ca which would be singular in an infinite system to the "normal" value of Ca for a degree 
of freedom with a correlation length I ~ 1 fm as 



(^)max ( £max\ V 



(Ca) norm \ » 



36 . (113) 



This does not mean that the specific heat Ca is multiplied by 36, because it receives a 
nonsingular contribution from other degrees of freedom. However, it suggests that in using 
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Eq. (|59D to estimate how much the anti-correlation induced reduction of Ft is weakened at 
the critical point, it is reasonable to expect that Ca can be up to an order of magnitude 
larger there than it is near T c far from the critical point. 

We now return to our discussion of the effects of the long wavelength sigma fluctuations 
on the fluctuations of the pions, encoded in the microscopic correlator (^). We derived ([88]) 
using mean field theory, and would now like to discuss the effect of non-mean-field corrections. 
We mentioned previously that fluctuations of the sigma field around the minimum of Q(a) 
are not small; we now argue that this does not make much difference to the quantity of 
interest. One way to see how these corrections can appear is to realize that, at higher order 
in G, diagrams with a-bubbles which are actually logarithmically divergent as m a — > 
will contribute. These bubbles have to be resummed and may modify the exponent of 
the m~ 2 singularity in ([8^). This exponent is easy to infer from universality arguments. 
Diagrammatically, the 1/m 2 is the zero momentum value of the sigma propagator, i.e., the 
sigma field susceptibility. For the 3d-Ising universality class we know the corresponding 
exponent to be 7 / v = 2 — rj which is ~ 2 to within a few percent because 77 is small. We can 
therefore safely use the mean-field formula flgSp with its m~ 2 divergence for our estimate, 
and will take m a ~ l/£ m ax ~ 1/(6 fm). It therefore turns out that even though the effects of 
Section 4 depend on Cy and the effects of this Section depend on the sigma susceptibility, 
both are controlled by the exponent 7/1/ « 2. 



5.5 Magnitude of the Effects 



We now have all the ingredients in place for our estimate of the size of the effect of the critical 
fluctuations of the sigma field on the fluctuations of the observed pions, via the coupling G. 
We reproduce here Eq. (|89D 



<(A<7) 2 > 




1 G 2 
Vm 2 T 



to 



(114) 



which we now apply for q = p?. We have restored the factor of V. The first term in the 
curly brackets includes the single particle inclusive distribution enhanced by the Bose effect. 
The second term is the effect we are interested in now. As we did in our estimate of the 



effects of energy conservation and thermal contact in Section [4.3| , we shall express the size 
of the effect of interest as the ratio of the entire expression in curly brackets in ( \1 14j ) to the 
first term in these brackets. We find 



1 + 0.35 



(Cf freeze— out ^ ( 'Cfrcezc— out 
300 MeV 
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for ^ = 60MeV 



for fi n = 



(115) 



(116) 
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where we have taken T = 120 MeV. As in Section 4, the effect will be diluted by about a 
factor of two because not all of the pions which are observed are direct. We have written the 
estimates ( |115| , |116| ) in such a way that the largest uncertainties are manifest. The size of 
the effect depends quadratically on the coupling G. We argued above that G is reduced to 
Ge ~ 300 MeV at the critical point but, as we explained, there are caveats in this argument. 
Furthermore, freeze-out may occur somewhat away from the critical point, in which case G 
would be somewhat larger, although still much smaller than its vacuum value. The size of 
the effect also depends quadratically on the sigma correlation length at freeze-out, and we 
have seen that there are many caveats in an estimate like £f ree ze-out ~ £max ~ 6 fm. Finally, 
the effect is sensitive to the value of p^. There are reasons to believe that p n may be smaller 
near the critical point than far from it. Recall that p^ is zero at chemical freeze-out, and 
then grows until thermal freeze-out. At the critical point, the transition temperature T c is 
somewhat lower than at lower baryon chemical potential p, and this suggests that T c h may be 
lower than the value measured in 160 AGeV collisions. Furthermore, we have argued in Ref. 
HJ that the thermal freeze-out temperature Tf will be somewhat higher in the vicinity of the 
critical point, because the system lingers there and expands for a while with a temperature 
near T c .0 If the temperature window T c h — Tf is small near the critical point, then p^ may 
be significantly smaller than 60 MeV there. 

We have studied two different effects on y/F in Sections 4 and 5. The effects of energy 
conservation and thermal contact lead to a few percent reduction in y/F ; this reduction will 
be significantly lessened near the critical point due to the divergence of Ca- This yields 
perhaps a 1 — 2% increase in y/F near the critical point. The direct effect of the critical 
fluctuations which we have estimated in this section is a further, larger, increase in y/F by 
a factor of y/F^. We have displayed the various uncertainties in the factors contributing 
to our estimates ( |115| , |116j ) so that when an experimental detection of an increase and then 



subsequent decrease in y/F occurs, as control parameters are varied and the critical point 
is approached and then passed, we will be able to use the measured magnitude of this 
nonmonotonic effect to constrain these uncertainties. It should already be clear that an 
effect as large as 10% in yfFZ is easily possible; this would be 50 times larger than the 
statistical error in the present data. 

Once we have evaluated the microscopic correlator (An p Auk) , we can estimate the effect 
of the critical fluctuations of the sigma field on the fluctuations of many different pion 
observables. To this point, we have focussed on ((Apt) 2 )', we now give a brief account of the 
effect on ((AN) 2 ) and (AN Ap?) ■ We can use the microscopic correlator ( |5SD to obtain 

,21 2 



((AN) 2 ) = 5> p 2 + 



1 G 
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for q = pr, T = 120 MeV and fi n = 0. The coefficient which is 1.0 in the last line of (|117 ) 



increases to 2.0 if p,^ = 60 MeV. We see that there can easily be a very large increase in 



21 r 



The expansion rate if of (111) does not decrease. However, the rate of change of T with time is reduced 



because of the large specific heat. 
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the multiplicity fluctuations of the direct pions near the critical point, due to the coupling 
between the direct pions and the critical fluctuations of the sigma field. As we have noted 
previously, the noncritical fluctuations of the total pion multiplicity are dominated by the 
pions from resonance decay. Using the result fl32|) , we estimate that the sigma induced 
critical multiplicity fluctuations of the direct pions can easily lead to a 10% increase in the 
total ((AN) 2 ). This is comparable in magnitude to the effect on ((Ap?) 2 ), and should be 
easily detectable. We will see in Section 6 that there should be even a further increase in 
the fluctuations of the multiplicity of those pions with low pr- 

Turning now to the cross correlation between an intensive quantity and N, we use the 
microscopic correlator fl88|) to calculate 

(ANAg) _J_ly nHq _ {q)) + J_9l\y<] \y v l iq _ {q)) 

(N)(q) it MV {q)) Vml T [\ u k \ [\ ^ {q)) 
/r \ 2 ( t ^ r 

I ^freeze— out \ / Sfrcczc— out 




-0.021 {1 + 12. nuuzc ~ our ^ ncczu ~ OUI ^ for fi n = 
* V 300 MeV J \ 6 fm ' ' 

« -0.046 |l + 13. (|=^) 2 ^Y^) 2 \ fOT ^ = 60 MeV ' (118 

where we have taken q = Pt and T = 120 MeV. This correlation only receives contributions 
from nontrivial effects, and we see that near the critical point, the contribution from the 
interaction with the sigma field is dominant. Once again, we see that this correlation is a very 
interesting quantity to use to look for the critical point. As the critical point is approached, 
thermal contact with a heat bath whose heat capacity is diverging reduces the effect of 
energy conservation, as seen in (|70D ; we now see that this reduction is overcompensated 
by the larger increase in the cross correlation which is induced by the direct coupling of 
the pions to the sigma field. The lesson is clear: although this correlation is small, it may 
increase in magnitude by a very large factor near the critical point. 

The effects of the critical fluctuations can be detected in a number of ways. First, one 
can find a nonmonotonic increase in F a , the suitably normalized increase in the variance 
of event-by-event fluctuations of the mean transverse momentum. Second, one can find a 
nonmonotonic increase in ((AN) 2 ). Both these effects can easily be between one and two 
orders of magnitude greater than the statistical errors in present data. Third, one can find 
a nonmonotonic increase in the magnitude of (AptAN) . This quantity is small, and it has 
not yet been demonstrated that it can be measured. However, it may change at the critical 
point by a large factor, and is therefore worth measuring. In addition to effects on these and 
many other observables, it is perhaps most distinctive to measure the microscopic correlator 
(An p Ank) itself. The term proportional to 1/m 2 in fl58| ) has a specific dependence on p and 
k. It introduces off-diagonal correlations in pk space. Like the off-diagonal anti-correlation 
discussed in Section 4, this makes it easy to distinguish from the Bose enhancement effect, 
which is diagonal in pk. Near the critical point, the off-diagonal anti-correlation vanishes 
and the off-diagonal correlation due to sigma exchange grows. Furthermore, the effect of o 
exchange is not restricted to identical pions, and should be visible as correlations between 
the fluctuations of tt + and 7r~. The dominant diagonal term proportional to 8 p k in (^) will 
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be absent in the correlator (An+ An k ), and the effects of a exchange will be the dominant 
contribution to this quantity near the critical point. 



6 Pions From Sigma Decay 

There is another signature of freeze-out near the critical point discussed in |IJ in addition 
to those we have analyzed in depth above. For choices of control parameters such that 
freeze-out occurs at or near the critical endpoint, the excitations of the sigma field, sigma 
(quasi)particles, are nearly massless at freeze-out and are therefore numerous. Because the 
pions are massive at the critical point, these cr's cannot immediately decay into two pions. 
Instead, they persist as the temperature and density of the system further decrease. During 
the expansion, the in-medium a mass rises towards its vacuum value and eventually exceeds 
the two pion threshold. Thereafter, the cr's decay, yielding a population of pions which do 
not get a chance to thermalize because they are produced after freeze-out. Here, we estimate 
the momentum spectrum of these pions produced by delayed a decay. An event-by-event 
analysis is not required in order to see these pions. The excess multiplicity at low will 
appear and then disappear in the single particle inclusive distribution as control parameters 
are varied such that the critical point is approached and then passed. 

The event-by-event fluctuations of the multiplicity of these pions reflect the fluctuations 
of the sigma field whence they came0. We estimate the resulting increase in the event-by- 
event fluctuations of N, the total pion multiplicity. 

We begin with the inclusive single particle pr-spectrum of the pions from sigma decay. 
We use the expression (|3|) for the width of the a, but now treat m a as time-dependent. 
We should also take G to evolve with time. However, the dominant time- dependent effect is 
the opening up of the phase space for the decay as m a increases with time and crosses the 
two-pion threshold. We will therefore treat G as a constant. In Section 5, we estimated that 
in vacuum with m a = 600 MeV, the coupling is G ~ 1900 MeV, whereas at the critical end 
point with m a = 0, the coupling is reduced, perhaps by as much as a factor of six or so. In 
this section, we need to estimate G at the time when m a is at or just above twice the pion 
mass. We will use G ~ 1000 MeV, recognizing that we may be off by as much as a factor of 
two. 

Let us parameterize the time dependence of the sigma mass by 

m a (t) = 2m n {l + t/r) (119) 

where we have defined t = to be the time at which m a has risen to 2m n and have introduced 
the timescale r over which m a increases from 2m n to Am^. We will be interested in times 
< t < r, for which this linear parameterization of the time dependence is not unreasonable. 
Note that with this choice of notation, freeze-out occurs at a negative time, and the collision 
begins at an even more negative time. Substituting into (f?3|), and working to lowest order 
in t/r, we find 

T ^~^^ = D ^< < 120 > 
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where we have defined 

D = j^- ~ (300 MeV) f t— - T ) . (121) 



32vrm^ v ' V 1000 MeV 

N(t), the number of cr's present at time t, is determined by 
1 dN(t) 



N(t) dt 
and is therefore 



-r(t) = -DJt/r , (122) 



iV(t) = JV(0) exp D t 3 / 2 r- 1/2 ) . (123) 

We can now estimate the momentum distribution of the pions produced in the decay 
of the sigmas, upon making the assumption that the sigmas are at rest when they decay. 
This is a good approximation for two reasons. First, as the system expands after freeze-out, 
the sigma mass is increasing as we have discussed. This means that the kinetic energy of 
each sigma is decreasing. Second, during the time between freeze-out and decay, some of the 
sigmas which happen to be moving outwards toward the less dense region of the collision 
in which their mass would increase more than allowed by energy conservation will instead 
be reflected back inward. Each sigma which suffers such a reflection loses momentum, as 
the reflection occurs as if off an outward moving surface. This effect confines the sigmas to 
the densest region of the plasma, where their mass remains low for the longest time, and in 
addition reduces their momenta. We do not attempt a quantitative estimate of these two 
momentum-reducing effects here. Suffice to say that since at freeze-out the typical sigmas 
will have momenta of order the freeze-out temperature or less, we think it reasonable to 
approximate them as being at rest at time t — when they begin to decay. 



Sigmas which decay at rest at time t each yield two pions with momenta p ~ m vr y2t/r, 
to lowest order in t/r. As a result, the number of pions with momenta m^Jltjr < p < 



m^y 2(t + dt)/r is —2dt(dN(t)/dt) with N(t) given by ( |123| ). Upon making suitable substi- 
tutions, we find that the number of pions with momenta between p and dp is 

dN = V2N(0)Dr P ^-ex P ( ^Dr^r) . (124) 

V ; ml \ 3^2 ml) K 1 

With the momentum distribution in hand, we determine the mean pion momentum to be 

2V6 3 -2/3 r ( 1 / 3 ) mx (D T )-i/3 = 1.45 m^CDr)- 1 / 3 . (125) 



p mc 



Large r corresponds to slow expansion and a sigma mass which consequently increases only 
slowly with time; large D corresponds to a large coupling constant G. It therefore makes 
sense that if Dt is large, the sigmas decay before the sigma mass has increased far above 
threshold, and the resulting pions have small momenta. We defined r to be the time it takes 
the m a to increase from 2m n to 4m n . This timescale is hard to estimate, but our result is 
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not strongly dependent on r. It seems likely that 5 fm < r < 20 fm and we therefore quote 
our result as 

_ inc n ^ Q /1000 MeV\ 2 / 3 /l0fm\ 1/3 

p inc ~0.58m^ J I — —J , (126) 



where we have used ( p.21| ). We therefore estimate that if freeze-out occurs near the critical 
point, there will be a nonthermal population of pions with transverse momenta of order half 
the pion mass distributed according to ( |124j) . 

How many such pions can we expect? That is, how large is iV(0)? This is determined by 
the sigma mass at freeze-out. If m a is comparable to at freeze-out, then there are half as 
many cr's at freeze-out as there are charged pions. Since each sigma decays into two pions, 
and two thirds of those pions are charged, the result is that the number of charged pions 
produced by sigma decays after freeze-out is 2/3 of the number of charged pions produced 
directly by the freeze-out of the thermal pion gas. Of course, if freeze-out occurs closer to 
the critical point at which m a can be as small as (6 fm) -1 , there would be even more sigmas. 
We therefore suggest that as experimenters vary the collision energy, one way they can 
discover the critical point is to see the appearance and then disappearance of a population 
of pions with (pt) ~ rn n /2 which are almost as numerous as the direct pions. Yet again, it 
is the nonmonotonicity of this signature as a function of control parameters which makes it 
distinctive. 

As we discussed briefly in |[J , the event-by-event fluctuations in the multiplicity of these 
low momentum pions are also of interest. If we were able to measure the multiplicity of 
sigma quasiparticles at freezeout, we would find fluctuations given by 

<(AAH 2 > = £K> (l + K>) , (127) 

p 

where the r^'s are the sigma occupation numbers. In the present analysis, we neglect the 
effects of interactions among the sigmas and just take as for an ideal Bose gas with small 
m a . We expect that this makes our prediction for the fluctuations an underestimate. Since 
m a is small, the low momentum modes have large occupation number, and have fluctuations 
proportional to the square of their occupation number. Each sigma eventually decays into 
two pions, whose momenta are determined by the time at which the sigma decays, rather than 
by the sigma momentum at freeze-out. It is therefore not possible to make a measurement 
on the pions which restricts the J2 P i n ( P-27|) to low pt- We therefore do the entire sum, and 
find that the variance of the event-by-event distribution of the multiplicity of the cx-produced 
pions is 



((ANY) = 2(N) (I + (njpj (128) 

where N is the number of charged pions. The factor of two arises because every sigma which 
produces charged pions produces two charged pions, and was discussed in Section 3. Taking 

— inc 

m a = yields (rip) ~ 0.37, and therefore 

((AN) 2 ) « 2.74(iV) . (129) 
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We have already seen in Section 5 that the critical fluctuations of the sigma field increase 
the fluctuations in the multiplicity of the direct pions sufficiently that the increase in the 
fluctuation of the multiplicity of all the pions will be increased by about 10%. We now 
see that in the vicinity of the critical point, there will be a further nonmonotonic rise in 
the fluctuations of the multiplicity of the population of pions with (pr) ~ ^?r/2 which are 
produced in sigma decay. 

7 Summary and Outlook 

In order to estimate the magnitude of the effects of critical fluctuations, one must first 
analyze the background, noncritical fluctuations. NA49 data from PbPb collisions at 160 
AGeV shows that the event-by-event distribution of the mean transverse momentum is as 
perfect a Gaussian as the central limit theorem allows. Since a system in thermodynamic 
equilibrium exhibits Gaussian fluctuations, in Section 3 we give a quantitative answer to the 
question of how much of the observed fluctuations are thermodynamic in origin. To this 
end, we model the matter at freeze-out as an ideal gas of pions and resonances in thermal 
equilibrium, estimate the resulting fluctuations, and compare with the data. 

We calculate the event-by-event fluctuations of px, an intensive quantity which is there- 
fore little affected by nonthermodynamic fluctuations in the initial size of the system. We 
find that the resonances turn out to be of little importance — the resonance gas prediction 
for ((Apt) 2 ) / (pr) 2 is almost indistinguishable from that of an ideal Bose gas of pions. Fur- 
thermore, we have verified quantitatively that the correlations between pions introduced by 
the fact that some originate in resonance decays can be neglected. We have computed the 
effects of Bose enhancement, and find that they increase ((Ap T ) 2 ) by only a few percent, 
although the precision of the data should make effects of this magnitude detectable. 

The difficulty comes in the treatment of the collective flow. This hydrodynamic expansion 
boosts the momenta of the pions, affecting both the numerator and the denominator in 
((Ap T ) 2 )/(p T ) 2 . Although we do expect that the effect cancels in the ratio to a significant 
extent, the "blue shift" approximation which we have used is too simple. We have shown 
quantitatively that the fluctuations in the flow velocity can be neglected. However, the effects 
of the flow itself are not sufficiently accurately treated as a uniform blue shift, and must be 
treated more quantitatively in the future. We find that our prediction for ((Ap T ) 2 ) / (pr) 2 is 
about 90% of that which NA49 observes. This gives us further confidence that we can use 
thermodynamics to understand the great bulk of the observed fluctuations; improving the 
precision of the prediction by improving upon the blue shift approximation remains to be 
done. 

The data are precise enough that we can do more than analyze the "bulk" of the fluc- 
tuations. We can ask, for example, about the ratio y/~F of (-/V) 1/,2 i> ebe (pT) to the variance 
of the inclusive single particle distribution. This ratio is insensitive to the effects of the 
flow velocity. For a classical ideal gas, \f~F~ = 1. We estimate that Bose effects result in 
y/F = « 1-02. In the data, however, VF = 1.002 ± 0.002. The Bose effects may be 

small, but they are ten times larger than the statistical error in the data. The Bose corre- 
lations are being compensated by some anti-correlation, and in Section 4 we find a possible 
explanation. 
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We show that energy conservation results in an anti-correlation which is reduced by 
thermal contact between the direct pions and an unobserved heat bath. The anti-correlation 
vanishes if the heat bath has infinite heat capacity. This effect, and indeed everything 
about the fluctuations we analyze, can be derived from the correlator {An p Arik) between 
the fluctuations of the occupation numbers of pion modes with momenta p and k. Energy 
conservation implies that if n p fluctuates up, then n k is more likely to fluctuate down. The 
magnitude of the effect depends on the heat capacity of the "heat bath", but we estimate 
that it leads to \/~F = \[Y^F^ with ^fF^ « 0.99. 




With more detailed experimental study, either now at the SPS, or soon at RHIC (STAR 
will study event-by-event fluctuations in p T) N, particle ratios, etc; PHENIX and PHOBOS 
in N only) it should be possible to disentangle the different effects we describe. Making 
a cut to look at only low px pions should increase the effects of Bose enhancement. Bose 
enhancement effects are sensitive to fi n , and measuring these effects would allow one to 
make an experimental determination of this quantity. The anti-correlation introduced by 
energy conservation and thermal contact is due to terms in (An p Arik) which are off-diagonal 
in pk. Thus, a direct measurement of (An p Arik) would make it easy to separate this anti- 
correlation from other effects. The cross correlation {AN Apr) is a very interesting observable 
to study because it only receives contributions from interesting effects, like Bose enhance- 
ment, thermal contact and the critical fluctuations discussed in Section 5. We hope that the 
combination of the theoretical tools we have provided and the present NA49 data provide a 
solid foundation for the future study of the thermodynamics of the hadronic matter present 
at freeze-out in heavy ion collisions. 

We also consider fluctuations in multiplicity N, an extensive observable. These are not 
affected by the boost which the pion momenta receive from the collective flow, and this makes 
them easier to calculate than the fluctuations in p T . However, multiplicity fluctuations are 
contaminated experimentally by fluctuations in the initial state, for example due to the 
distribution of impact parameters. This experimental contamination can be reduced by 
making a tight enough centrality cut using a zero degree calorimeter. We compare the 
multiplicity fluctuations of the 5% most central events in the NA49 data to those we predict 
from a resonance gas, and find evidence that about 75% of the observed fluctuation is indeed 
thermodynamic in origin. We find that resonances play a significant role in this comparison, 
increasing the magnitude of thermodynamic fluctuations of the pion multiplicity and bringing 
it closer to the data. 

With the foundations established, we then describe how the fluctuations we analyze will 
change if control parameters are varied in such a way that the baryon chemical potential at 
freeze-out, /if, moves toward and then past the critical point in the QCD phase diagram at 
which a line of first order transitions ends at a second order endpoint. We provide quantita- 
tive estimates of the magnitude of the change in the observables we have analyzed which can 
be expected near this point. The agreement between the noncritical thermodynamic fluc- 
tuations in pt which we analyze in Section 3 and NA49 data make it unlikely that central 
PbPb collisions at 160 AGeV freeze out near the critical point. Estimates we have made in 
a previous paper suggest that the critical point is located at a baryon chemical potential fi 
such that it will be found at an energy between 160 AGeV and AGS energies. This makes it 
a prime target for detailed study at the CERN SPS by comparing data taken at 40 AGeV, 
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160 AGeV, and in between. If the critical point is located at such a low /i that the maximum 
SPS energy is insufficient to reach it, it would then be in a regime accessible to study by the 
RHIC experiments. We want to stress that we are more confident in our ability to describe 
the properties of the critical point and thus to predict how to find it than we are in our 
ability to predict where it is. 

The critical fluctuations near the endpoint affect the event-by-event fluctuations which 
we analyze in two different ways. First, all effects of energy conservation should be greatly 
reduced by thermal contact as the critical fluctuations in the sigma field cause the heat 
capacity to grow. Second, these critical fluctuations have direct effects on the fluctuations of 
the pions through the Gaim coupling. We analyze the most singular effects of this coupling, 
which are due to the zero momentum mode of the sigma field. It is possible to analyze 
subleading corrections using a diagrammatic approach, but we leave this to the future. 

In the chiral limit, the critical point becomes a tricritical point at which G vanishes. We 
estimate the vacuum value of G and use scaling laws valid near a tricritical point to estimate 
G at the critical point. We then estimate the increase in the fluctuations of N and pr 
distributions which we expect in heavy ion collisions which freeze out near the critical point. 
Finite size and finite time effects prevent ((AN) 2 )/(N) and ((Apt) 2 )(N) from diverging, as 
they would in an infinite system. We estimate that ((AN) 2 )/(N) can grow by more than 
10%. The ratio VF which describes the pr-Auctuations becomes \J F B F T F a with \[F\, about 
1.1. This effect is not large but is still predicted to be 50 times larger than the statistical 
error in the present NA49 measurement of y~F = 1.002 ±0.002. We quantify the uncertainty 
in our estimates in terms of the sigma correlation length £ and the coupling G at freezeout; 
measurement of the enhanced fluctuations of N and pt would allow one to estimate G£. 

We want to emphasize that the ratio y/F is not the only observable which can be used 
to detect the proximity of the critical point, and indeed is not the most sensitive observable 
available. We have focussed on \fF because it is simple to define, and because NA49 has 
published data to which we can compare our predictions. However, the specific form of 
the singularity in (An p Arik) which we find in Eq. flS8|) tells us how to construct observables 
which are more sensitive to the critical fluctuations. One possibility is to consider observables 
which are sensitive to the off-diagonal part of (An p Ank) , because the noncritical off-diagonal 
anticorrelation in (An p Ank) should be replaced by a much larger off-diagonal correlation 
near the critical point. A second possibility is an analysis of the cross correlation (AN Apt)- 
Because this cross correlation is dominated by interesting effects, we have seen that it can 
increase by an order of magnitude at the critical point. A third possibility is to construct 
a ratio like yF, but using only soft pions, with px less than a specified cutoff. The effects 
of the critical fluctuations are largest on the softest pions, and they are therefore masked in 
y/F which receives significant contribution from harder particles. Whereas we have found 
that the critical fluctuations change y/~F at the 10% level, their effect on a "soft \/F" can 
easily be at the factor of two level. 

Although the sigma quasiparticles themselves cannot be reconstructed, their presence can 
be detected even more directly than via their influence on the pions at freeze-out. If freeze- 
out occurs near the critical point, some time after freeze-out the sigma mass rises above the 
two pion threshold, and the sigmas decay quickly. Since these pions do not rethermalize, the 
resulting excess in the low-pr region of the pion momentum spectrum should be observable. 
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The mean p? of these pions is about m 7r /2, and they are almost as numerous as the direct 
pions. The event-by-event fluctuations in the multiplicity of these soft pions would be even 
larger than those of the rest of the pions near the critical point. 

In summary, our understanding of the thermodynamics of QCD will be greatly enhanced 
by the detailed study of event-by-event fluctuations in heavy ion collisions. We have esti- 
mated the influence of a number of different physical effects, some special to the vicinity 
of the critical point but many not, on the fundamental correlator (AnpAn^) . This is itself 
measurable, but we have in addition used it to make predictions for the fluctuations of ob- 
servables which have been measured at present, like {(Apt) 2 ) and {(AN) 2 ) and also for the 
cross correlation (ANAp?)- The predictions of a simple resonance gas model, which does 
not include critical fluctuations, is to this point in broad agreement with the data. More 
detailed study, for example with varying cuts in addition to new observables, will help to 
further constrain the nonthermodynamic fluctuations, which are clearly small, and better un- 
derstand the different thermodynamic effects. The signatures we analyze allow experiments 
to map out distinctive features of the QCD phase diagram. The striking example which we 
have considered in detail is the effect of a second order critical end point. The nonmonotonic 
appearance and then disappearance of any one of the signatures of the critical fluctuations 
which we have described would be strong evidence for the critical point. Furthermore, if a 
nonmonotonic variation is seen in several of these observables, then the maxima in all the 
signatures must occur simultaneously, at the same value of the control parameters. Simulta- 
neous detection of the effects of the critical fluctuations on different observables would turn 
strong evidence into an unambiguous discovery. 
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Note added in proof 

As we have stressed in Section 3, the fluctuations in an extensive quantity such as the 
observed multiplicity are unlike fluctuations in intensive quantities in that they receive 
significant contributions from both (i) thermodynamic fluctuations at freeze-out and (ii) 
nonthermodynamic fluctuations during the initial stage of the collision. Our approach has 
been to use a comparison between the data and thermodynamic predictions to constrain the 
magnitude of non-thermodynamic fluctuations. 

After this paper was submitted, Ref. |33| appeared. These authors have attempted a 
theoretical treatment of those non-thermodynamic fluctuations which are purely geometrical 
in that they can be attributed to the distribution of impact parameters. Further analysis is 
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presented in Ref. [33|. These authors include in addition the effects of fluctuations in the NN 
cross-section [[H]], which they find to be small, and also fluctuations in the number of punch- 



through spectators and effects due to the diffuse edges of the incident nuclei, both of which 
are significant. Combining all contributions to the multiplicity fluctuations, thermodynamic 
and nonthermodynamic, yields fluctuations which, with no new free parameters, reproduce 
the magnitude of the observed multiplicity fluctuations to within a few percent accuracy 
|31|. 



A Finiteness of Multiplicity 

Throughout the body of the paper, we use event-by-event and inclusive averages defined 
probabilistically If we were interested in an infinite ensemble in which each member of the 
ensemble was in the thermodynamic limit, no translation would be required. However, when 
we want to compare the relations involving quantities which are defined probabilistically 
to those measured in an experiment, as in any application of probability theory we must 
have estimators for these quantities which can be constructed from finite samples. In this 
Appendix, we discuss the effects due to finite sample size. 

The typical size of these effects is one over a power of the sample size. The total number, 
C, of events can be easily made very large (say, 10 6 ), so that 1/C, and even 1/yC is much 
smaller than the physical effects we consider, (such as Bose enhancement, for example) which 
are of the order of a few percent or more. However, the number of pions in an event, N, is 
limited by the size of the colliding system and the experimental acceptance of a detector, 
and is typically of the order of a few hundred. This can introduce corrections of the order 
of a fraction of a percent. Of course, these effects are negligible when compared to the bulk 
of the fluctuations, which we analyze in Section 3. They are also smaller than the effects 
we discuss in Sections 4 and 5, where we are interested in signatures which rise and fall by 
of order 10% near the critical point. However, we have seen in Section 3 that the statistical 
errors in the present data are small enough that one can compare quantities like {N)v^ he and 
v? QC to a precision of less than a percent. At this level, we must understand how to deal with 
the 1/N corrections. 

Let us consider a sample of values of some one-particle observable q. This sample is 
broken into C subsets, i.e., events, with N a values per event. We use the notation: qf, where 
a = 1 . . . C and i = 1 . . . N a . (For example, qf may be the momentum of the z'th pion 
in the a'th event.) The numbers qf are distributed according to some (joint) probability 
distribution. We assume that the expectation value is the same for all qf: M[qf] = m. 

In Section 2 we pointed out that the inclusive mean q mc is the same as the event-by-event 
average (q). Both quantities are defined in the infinite C and N limit. Let us now try to 
estimate these two quantities, using our finite sample qf. The natural estimate for q mc is the 
following: 



V C v iV a a 1 C N a 

Z ^q=l Z-a=l Hi 1 \ -> \ ~ 



Ct = ^ a =c : - =7™LL?i ( 13 °) 
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where we have introduced (somewhat inconsistently, but suggestively): 



(A0 = ^X> a , (131) 

° 0=1 

which is (an estimate for) the mean multiplicity in an event. The total number of q's in the 
sample is C(N). The property of the estimate ( |130| ) is that its expectation value is equal to 
m: 

M[q™} = m, (132) 



for any N or C. The standard deviation of this quantity is 0(l/yC(N)), so the estimate 
becomes perfect in the infinite C, N limit. Now let us estimate the event-by-event mean of 
q. An estimate which appears natural is: 

(,> -= csUg'O- (133) 

The expectation value of this estimate is: 

M[(g) cst0 ] = m, (134) 
and this estimate also becomes perfect as C, N go to infinity. However, 

(?U^. (135) 

One can show that the difference between the two is on the order of ((AN) 2 ) /N 2 ~ 1/N. It 
is obvious how to improve the estimate (|133| ) to make the relationship (q) = q mc hold exactly 
for finite N. Writing: 

1 c N ( 1 Na \ 

Ct c = 7; E 7^ iw E = to)*> ( 136 ) 



we can interpret this definition of the estimate (g)est as a result of averaging over events with 
a weight proportional to the multiplicity in this event, N a . It is also clear intuitively that 
such a procedure is more natural than taking all events with equal weight as is done in ( |133| ) . 
What is important is that, by construction, this procedure rids us of any 1/N correction to 
the equality between (q) est and 

Let us now consider estimating variances of q. The natural estimate for the inclusive 
square variance is: 

1 c N a 

U \ iV / _ 1 0=1 i=l 

Assuming that the variables g" are uncorrelated and their dispersions are equal, i.e.: 

M[qU b i} =m 2 + a 2 5 tJ 5 a \ (138) 
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one can show that 



M[vUlU\ = °\ (139) 

for any C, or N. This is the property which we require of this estimate. 

Next, we consider estimating the event-by-event variance. One seemingly natural candi- 
date is: 

1 c ( 1 Na \ 2 

vL(qUo = E [w E - (iUJ • ( 14 °) 

This is the procedure used by NA49 to calculate f e beG°T) from their data, leading to the 
result shown in Table |2|. Let us calculate the expectation value of this quantity, assuming 
that all g's are independent as in (|138|). We find: 



MKJ,U„) = ^ (^) - ^) « (141) 

where in the last approximate equality we neglected terms of relative size 0(1/C). In the 
case of completely uncorrected q we expect the following relation to hold between the v^ he (q) 
and vf nc {q): 

(N)v 2 eh M=vlM- (142) 

This equality cannot and should not be satisfied for an arbitrary sample (unlike the equality 
(q) = q mc ), but we want it to be satisfied on the level of expectation values: 

(A0M[t; e 2 be (g)csto] = M[vUqU]. (143) 

We see that for the estimate (|140|) the difference between the l.h.s. and the r.h.s of eq.( |143| ) 
is 

(N)M[vl e (qUo] - M[vlMest] « a (V>(^) - l) « ^ 2 ^P, (144) 
where we have neglected the 0(1/ C) corrections and the corrections which are higher order in 



1/N. The difference ( |144j ) is a 1/N effect, of course, but with iV ~ 300 it could easily reach a 
fraction of a percent. Note that we are not talking here about statistical fluctuations around 
the mean values which make the two estimates deviate from sample to sample. Such effects 
are of the order 1 / \J~C and can easily be made negligible with sufficient statistics. Equation 
( |144j) , on the other hand, reflects a systematic discrepancy between the expectation values 
of the estimates, which would only go away if N, the number of particles in one event, were 
infinite. 

Our task now is to give an estimator for f| be (g) which satisfies ( |143j ) with no errors of 
order 1/N in the case when the qf are uncorrelated. The lesson we learned from the estimator 
for (q) suggests that we take each event with the weight N a /(N). This gives: 

c AT / 1 N a 



c 



C N ( ^ N a y 

t£m ' (145) 
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Calculating the expectation value we find that 



a 2 



M[vi e (qU] = (146) 

exactly! This means that the estimate (|145|) satisfies our criterion ( |143|) exactly, to all orders 
in 1/N and 1/C. This is the estimate for v 2 he (q) that should be used to analyze experimental 
data. It introduces no 1/N or 1/C errors in the statement that the ratio F = 1 in the absence 
of correlation or interaction between the pions. 



We can avoid having to apply the formula ( |145| ) to the original experimental data set in 
order to recalculate the f e be(j>r) given in Table |^, which used the estimate ( |140| ). Using our 
analysis, we can instead just use the fact that (the expectation values of) the two estimates 
are related by (see ( |144| )): 



, N ( 1 (( AN ) 2 )\ , X 

Webe(<?)cst = V chc {q) cst0 II - ~ ^ 2 I , (147) 

up to corrections which are higher order in 1/N (and corrections of order 1/y/C). We use 
this relation to convert from one estimate to another in Eq. Q4"2] ) of Section 3.4. 

To finish the discussion of the 1/N effects we also point out that yet another estimate 



for Vebe(q) is used implicitly in the definition of $ PT in p|, |21 



$ pt = (N) 1/2 v ehe (p T ) est ^ - v inc (p T ). (148) 

The definition of f e be(pr)est* corresponds, in the language that we use here, to giving each 
event a weight N 2 / (N) 2 (and using C instead of C — 1): 



'ebe( ( ?)est* 



1 c N 2 ( 1 Na \ 2 

cSwUg#-«-f (i49) 



Calculating the expectation value of this estimate, one finds: 

M[vL( g U*] =- 2 {j^~ 75^3) - ^Ty (150) 
where we neglected 0(1/ C) corrections in the last step. We see that in the C — ■> 00 limit this 



estimate does not suffer from 1/N corrections as far as the relation (|143|) is concerned, and 



does not differ from the estimate (|145 ) . However, it does introduce small 1/C corrections, 



while the estimate (|145|) satisfies ( |143j ) exactly. 
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